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SYNOPSIS 

HYPERSONIC RAREFIED FLOW IN THE 
STAGNATION REGION OF A BIUNT BODY 

by 

Vipparthi Adimurthy 

Submitted in Partial Fulfilment 
of the Requirements for 

Ph.D. Degree 
to the 

Department of Mathematics 
Indian Institute of Technology, Kanpur 

July, 1972 

Wh en a blunt-nosed body descends the atmosphere at 
hypersonic speeds, an entire spectrum of flow patterns ranging 
from the free-molecule regime to the boundary-layer limit, 
passes around its stagnation point. So far, from the continuum 
point of view, the transition zone between these two limits 
is investigated either within the higher-order boundary-layer 
theory or with some thin-layer approximation to the Navier- 
Stokes equations. Both of these approaches have their limita- 
tions and cannot adequately describe the flow in the low-Reynol; 
number regime. In the present investigation, the full Navier- 
Stokes equations are used to describe the hypersonic rarefied 
gas flow near the stagnation point of a blunt body. Emphasis 
is placed to investigate the effect, of slip and temperature j 
jump boundary conditions on the structure and overall charac- j 
teristics of the flow. Using the concept of local similarity, j 
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the equations are reduced to an eighth-order system of non- 
linear, coupled ordinary differential equations with two point 
boundary conditions. For monatomic-and diatomic-gas streams, 
these equations are numerically solved with reasonable accuracy, 
with and without slip and temperature jump boundary conditions 
for adiabatic- and cold-wall cases. In the numerical computa- 
tions for different Mach numbers, the Reynolds number is varied 
in close steps over a wide range stretching from the fully- 
merged-layer to the boundary-layer regimes. 

It is found that slip and temperature jump boundary 
conditions at the surface affect the entire structure of the 
flow field (including the shock-wave-like region). These effects 
increase quantitatively as the Reynolds number is decreased. 

For adiabatic -wall case, the slip velocity increases to as 
much as 60/. of the corresponding free-stream velocity. Slip 
velocity and temperature jump have opposite effects on the flow. 
The total effect depends on the temperature of the body. Comp- 
ared to a cold wall case, slip effects are more on an adiabatic- 
wall. Also, effects of slip are more for a monatomic gas than 
for a diatomic gas. 

With a view to find the limit of validity of the thin- j 
layer approximation, a reduced form of the Ravier-Stokes 
equations also is integrated with and without slip and tempera- 
ture jump boundary conditions. Further, the present results are 
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oompared with the higher-order hound ary- layer theories and 
two-thin- layer-model of flows as well as with several experi- 
mental results. 

It is found that the thin-layer and two— layer-model 
theories fail to describe correctly the detailed flow field 
at low Reynolds numbers, even though they are able to predict 
surface characteristics like heat-transfer rate to the body in 
good agreement with experiments. Detailed flow comparisons 
with present full Navier-Stokes solutions show that the soluti- 
ons using the thin layer assumption are not valid below Re^lSO. 
Higher-order boundary-layer theories are found inadequate to 
describe the flow in such a low- Reynolds -number regime. 

Extensive comparisons of the present results with the ; 

merged-layer experimental density profiles (Russel, D.A, , 

Physics of Plaids, 1968) and temperature profiles (Ahouse, D.R, j 
and Bogdonoff, S.M. , AIM Pluid and Plasma Dynamics Conference, I 
1969) show a good agreement. Purther, the present results show 
a good agreement with the low-Reynolds-number impact-pressure 
measurements (Potter, J.L. and Bailey, A.B* , AIM Journal, 1964 ) I 
and recovery-temperature measurements (Hickman, R.S. and Giedt , : 

W.H. , AIAA Journal, 1963). The thin-layer theories fail to ; 

- . " ‘ . i ■ " * * , ~ i ' . ' ■ , ■: • ■ • ■ 

predict several of these experimental features even qualitati- ! 

' f 

vely. These factors establish the superiority of the full Navier- 
Stokes equations over the thin-layer theories in describing the 
characteristics of merged-layer flows. 
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CHAPTER 1 
INTRODUCTION 

1.1 Introduction 

The problem of predicting the aerodynamic character- 
istics of bodies flying at high speeds and high altitudes 
is of great importance in designing reentry vehicles. The 
flow patterns existing near the stagnation point of a blunt- 
nosed space-craft depend upon the degree of rarefaction in 
the atmosphere. As the body descends the atmosphere, the 

surrounding fluid changes its characteristics from the. free- 
molecule regime in the upper atmosphere to the continuum 
regime in the lower atmosphere. The free-molecule flow is 
well investigated from kinetic theory considerations. The 
other extreme of the hi^ily dense air, where the shock Wave 
is extremely thin and the viscous effects are confined to a 
narrow region near the body surface, lends itself to the 
classical, boundary-layer treatment . Between these two limits , 
we have situations in which the gas is sufficiently rarefied 
so that the appropriate mean free path becomes too large for 
the use of classical boundary-layer theory, but not large 
enough for the frce-molecule concepts. to apply. This .region 
may be called the intermediate regime. The main features of 
this regime are the thickening of the boundary layer and 



2 


the shock wave and their subsequent merging * This regime 
is usually subdivided according to the different theoretical 
models which are necessitated to describe the flow. The cha- 
racteristics of these regimes are discussed in detail in a 
later section. 

In general, a complete calculation of the flow field 
in the intermediate regime requires a full formulation of 
the problem based upon the kinetic theory approach. It had 
been considered until late fifties that some modifications 
of the Navier-Stokes equations, such as the Burnett equations 
or the Thirteen Moment equations, with slip velocity and 
temperature jump boundary conditions are applicable to this 
regime* But later experimental and theoretical findings for 
low-Mach-number shock-wave structure have indicated that the 
Navier-Stokes equations axe superior to the Burnett and Thir- 
teen Moment equations (e.g, . Sherman, 1955)., With, this view, 
in the last ten years a number of theoretical investigations 
of the rarefied-blunt -body problem are made using the Navier- 
Stokes approach. Note-worthy of these continuum. theories is 
the two-thin-layer theory of Cheng (1961, 1963a, 1966), which 
predicts the heat-transfer rates in a good agreement with . 
experimental results even at highly rarefied flow conditions* 
However, more recent experimental studies of the entire flow 
field (e.g* Russel 1968, Ahouse and Bogdonoff 1969) show$ 
that the existing continuum theories fail to predict the 
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detailed flow profiles* In our view, "this poor agreement in 
the detailed flow structure does not prove that the Navier- 
Stokes model is invalid in these situations. Probably, this 
failure occurs due to incorporating various approximations 
in the Navier-Stokes equations. In Section 1,2, we have 
discussed in a greater detail tho validity of using the 
Navier-Stokes equations to describe rarefied flows. The slip 
and temperature jump boundary conditions also are discussed. 

The existing theoretical investigations of the rare- 
ficd-blunt-body problem, from the continuum point of view, 
fall into three main categories. Firstly, we have the higher- 
order boundary-layer theories, developed systematically from 
the theory of singular perturbations. The classical boundary- 
layer theory of Prandtl is embedded in this theory as the 
first approximation. Higher-order equations offer corrections 
to this first approximation. So far, the second and third 
approximations arc worked out. The hi^ier-order boundary- 
layer theory is valid to describe the flow only, when there . 
is a small departure from the classical boundary- layer limit. 
The second approach to the blunt-body problem is the two- 
thin-layer theory, in which equations similar to Prandtl*s. 
boundary- layer equations are used to describe the flow betw^~ 
the shock wave and the body. Transport effects just behind, 
the shock wave are taken into account and consequently the 
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classical Rankine-Hugoniot shock conditions are modified. 
Investigations of this nature are found to he very successful 
in predicting surface characteristics such as the heat- 
transfer rate even upto the free-molecule limit, hut fail 
to describe the detailed flow field. The last and the most 
reliable approach to the problem is direct numerical integ- 
ration of the governing equations from the surface to the 
free stream. Unfortunately a thorough investigation of the 
rarefied-blunt -body problem using this approach is lacking 
possibly clue to the complexity of the problem and the comput-i 
tine involved in the numerical computations. A detailed, 
description of these three approaches, and their relative 
merits and demerits is presented in Section 1.4. 

The test of any meaningful theory lies in its being 
able to predict the characteristics of the. flow in agree-. . 
ment with experimental data. Unfortunately, until recently, 
the merged-layer regime of the rarefied gas dynamics remain- 
ed relatively immune to experimental flow-field investiga- 
tion. With recent development of sophisticated experimental 
techniques* like the electron-beam fluorescence and X-ray 
methods, an entirely new approach to the problem has become 
available, allowing detailed probing of the flow field. The 
available experimental results for the rar ef i ed-blunt -body 
problem are briefly discussed in Section 1.5* 
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The main aim of the present dissertation is to 
investigate thoroughly and extensively the rarefied, conti- 
nuum flow near the stagnation region of a blunt body under 
hypersonic conditions. By continuum flow it is meant that 
the Navier-Stokes equations are assumed to describe the flow. 
This investigation is motivated by the fact that the existing 
continuum theories are inadequate and limited in several 
respects. In this work, special emphasis is laid to under- 
stand the phenomenon of merging of the shock wave and the 
boundary layer, to determine the effect of slip velocity and 
temperature jump boundary conditions, and to compare our 
results as far as possible with available experimental results. 
In order to achieve these objectives, extensive numerical 
integrations are carried out. The aim and scope of the present 
investigation are detailed in Section 1.6. 

1 . 2 Discus sion on th e Na vier-Stokes Equations 
and the Slip Boundary Conditions 

The present investigation is primarily aimed at 
understanding the merged-layer-regime, in which the shock 
wave and the viscous layer near the body are thick and merge 
into each other. The formulation of valid differential 
equations in the merged-layer-regime is a matter still not 
resolved completely. A full kinetic-theory formulation for 
this regime is of prohibitive complexity. In such a formulation 
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one has to consider multiple collisions of the molecules* 

Further, the presence of a body in the flow highly compli- 
cates the analysis. It had been considered for some time 
that some modification of the Navier-Stokes equations, such 
as the Burnett equations or the Thirteen Moment equations, 
together with slip velocity and temperature jump boundary 
conditions could be regarded as applicable. Subsequently, 
this view had to be modified in the light of experimental 
as well as theoretical evidence. For example, we have the 
initial experimental studies of the shock structure caf” ' 
out by Sherman (1955), Sherman has studied the shock struc- 
ture at “1.8 in helium and l^ a 1,8 and 3,7 in air. He 
has found that at Ik - 1.8 the best fit to experiments is 
afforded by the Navier-Stokes theory with zero bulk viscosity 
in case of helium, and the bulk viscosity equal to 2/3 times 
the shear viscosity in case of air* At this Mach number, the 
profile obtained using the Burnett equations is very much 
different from that obtained using the Navier-Stokes equa- 
tions. It may be noted from Sherman's work that even taking aze-l 
1 --^nosity gives better agreement of the Navier- 
Stokes profile for air with the experimental profile as 
compared to the solution of Burnett equations* At Ik « 1 .8, 
the solution of the Thirteen Moment equations is not available, j 
but a tentative comparison of the available solution at I 

Ik « 1.6 shows that the Thirteen Moment equations fail to 
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predict the correct shock structure, giving a thicker shock 
than what is experimentally observed. But none of the the- 
ories, including the Navier-Stokes theory (with hulk visco- 
sity) affords good agreement with the shock structure in air 
at = 3.7. It is not clear how much importance one should 
attach to Sherman's result that the bulk viscosity is about 
two-thirds of the shear viscosity. The reasons for this are; 

1 ) this relation is inferred from the low -Mach-numb er shock 
structure comparisons, where even a zero bulk viscosity gives 
a better agreement of the Navier-Stokes profile with experi- 
mental shock structure, compared to the kinetic theories, 

2) experimental measurements of the bulk viscosity for li- 
quids (e.g. Roserihead, , 1954) show no correlation between the 
bulk viscosity and shear viscosity, and 3) the way Sherman 
defines the bulk viscosity it follows that bulk viscosity 
being equal to two-thirds of the shear viscosity would mean 
that the shear stress does not at all depend on the divergence 
of the velocity (from which after all the bulk property 
stems). After Sherman's work, Russell (1965), Camac (1965) 
and Robben and Talbot (1966) have measured the shock thick- 
nesses in high-Mach-number shocks, the last going upto 

^ = 17. Ih^y all find that the shock thicknesses predicted 
by the Navier-Stokes theory for high Mach numbers are much 
smaller than the experimental prediction. Russel (1965) finds 
by detailed comparison of the density profiles that on the 



8 


high-density side of the shock the discrepancy between the 
experiments and the Navier-Stokes theory is mot very much, 
though it becomes appreciable on the low-density side of, the 
shock. However, Russel (1965) and Camac.(1965) point out the 
importance of the dependence of the theoretical solution on 
the viscosity-temperature relationship, A correct viscosity 
law is found to bring the Navier-Stokes results much nearer 
to the experimental results. The Thirteen Moment and the 
Burnett equations are confined to monatomic gases. Further, 
some exact solutions of the Boltzmann equations for shearing 
flows (Ikenbury and Trusdell, 1956, Trusdell, 1956) show 
that the distribution function cannot be adequately approxi- 
mated by either the Burnett or the Thirteen Moment represen- 
tations. In these cases, the Navier-Stokes results are found 
to give better agreement with the exact solutions than the 
Burnett and Thirteen Moment results. 

With the above mentioned limitations of the Burnett 
and Thirteen Moment equations in view, in the past ten years 
a number of theoretical investigations of the rarefied blunt - 
body problem have been made using the Navier-Stokes model, 
as suggested by Adams and Probstein (1958)* Notable among 
these treatments are those of Cheng (1961, 1963a, 1966), Van 
Hyke (1962, 1963, 1964)i levinsky and Yoshihara (1962). and 
several others* All these approaches are critically discussed 
in Section 1.4, While each employes a different approach to 
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the problem, almost all of them use a thin-layer assumption, 
which means that the upstream disturbance created by the 
body is small in comparison to the nose radius. Recently, 
on the basis of comparison of experimental data with these 
approximate theories several investigators have strongly 
doubted the validity of the Navier-Stokes approach to the 
rarefied blunt-body problem, for example, Ahouse and Bogdon- 
off (1969) , who compare their experimental temperature 
profile with Cheng f s theory, and Vogenitz and Takata (1970), 
who compare their Monte Carlo profiles with the continuum 
profiles of Levinsky and Yoshihara (1962), come to such a 
conclusion. But these comparisons are not conclusive because 
the existing theories make several simplifying assumptions 
on the Navier-Stokes equations. It is one of the main purpo- 
ses of the present dissertation to investigate if the use 
of the full Navier-Stokes equations, without any thin-layer 
assumption, will make any difference to the situation. Indeed, 
it happens to make a significant difference and the full 
Navier-Stokes equations are found to predict the flow field 
in the fully-merged -layer regime with reasonable accuracy* 

In the merged-layer regime slip velocity and temper- 
ature jump boundary conditions are important rarefaction 
effects. The gas immediately adjacent to a solid surface 
possesses a finite velocity with respect to the surface. Also, 
the temperature of the gas next to a surface differs from the 
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surface temperature. These boundary condition effects become 
more important as the gas becomes more and more rarefied. 

For this reason, the merged-layer regime can also be called 
as the slip-flow regime. 

Several authors have derived the expressions for 
slip velocity and temperature jump boundary conditions using 
some approximate solutions of the Maxwell-Bo It zmann equation 
in the vicinity of the wall. The results given by Kennard 
(1938) are: 


u ■* 
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where x and y are respectively the coordinates along and 
normal to the surface of the body, u the tangential velocity* 
and Xyj the mean free path at the surface of the body, a is 
the thermal accommodation coefficient and cr the reflection 
coefficient. For the hypothetical cases of (i) entirely , 
specular reflection with vanishing energy exchange a = q =.0 
and ii) entirely diffuse reflection which is congpletely acco- 
mmodated to the surface temperature o= a « 1. .It may be. 
noted that even for entirely diffuse reflection and eomnleto 
thermal accommodation the slip velocity and temperature jump 
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do not vanish. The tern in the expression for 

clip velocity is known as the r theraal creep' term. According 
to this term a temperature gradient along a surface induces 
a flow in the direction of increasing temperature. This term 
vanishes when the temperature of the surface is kept constant, 
and can be neglected if the temperature variation along the 
surface is small in the region of interest. In general the 

y’ 

crebp term belongs to a higher order of smallness than the 
term containing (Maslen, 1952). As such, this term is usu- 
ally neglected in the studies of slip-velocity effect. The 
above expressions for the slip velocity and temperature jump 
envisage a small departure from equilibrium, that is, the 
velocity and temperature gradients should be small. Nonlinear 
corrections for large values of velocity and temperature gra- 
dients and for very low densities are suggested by Schamberg 
(1947) using the Burnett equations, and by Grad (1949) using 
the Thirteen Moment equations. But the usefulness of these 
equations itself is subject to doubt. Further limitations 
of the slip velocity and temperature jump boundary conditions 
are discussed by Schaaf (1963). 

In the present investigation we use the full Navier- 
Stokes equations with slip velocity and temperature jump 
boundary conditions (neglecting the creep term) to describe j 
the merged-layer flow in the stagnation region of a blunt j 

body. Unit reflection and thermal accommodation coefficients I 

/ : ; ; . ! 
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are used (a = a = 1). In view of the available experimental 
values of a and a (e.g,. Schaaf and Chambre', 1961), the 
present values of o and a appear to be representative in 
finding the effect slip velocity and temperature ;jump bou- 
ndary conditions on the flow field. In what follows, these 
boundary conditions are usually referred to simply as the 
slip boundary conditions* 

1 .3 Rarefied-Plow Regimes 

When a space-craft descends the atmosphere, an entire 
spectrum of flows ranging from the free-molecule flow to 
the boundary- layer limit passes around the stagnation point. 
At high altitudes, the air is rarefied and the motion of the 
individual gas particles is important. The parameter which 
determines the degree of rarefaction is the Enudsen number 
Kh, which is the ratio of an appropriate mean free path to 
a characteristic length in the flow. The delineation of the 
rarefied -flow regimes using the Ehudsen number concept is 
first given in a pioneering paper by Tsien ( 1946 ), where he 
suggests that either a body dimension like the nose radius 
or the boundary-layer thickness can be taken for the charac- 
teristic length. However, the use of the boundary-layer thic- 
kness, which comes from a theory valid for high Reynolds 
numbers, is not proper while dealing with rarefied flows, 
because the viscous flow may not be of the boundary-layer 
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type. The choice of an appropriate mean free path also has 
demanded considerable attention. A naive choice is to take 
\o. "k^-e mean free path corresponding to the free stream. 

But, as noted by Adams and Probstein (1958), the free-stream 
mean free path is not a characteristic mean free path of the 
flow. The density in the shock layer is appreciably more than 
that in the free stream, so that the mean free path near the 

surface of the body is much smaller than in the free stream. 
Probstein (1961) suggests that the appropriate mean free path, 

\ 0 , is that of the molecules emitted from the surface relative 
to the incident molecules. This mean free path is the smallest 
of the mean free paths that can be calculated among various 
combinations of the emitted and incident molecules. Assuming 
a very cold body and that the emitted molecules have the same 
temperature as the body itself, Probstein (1961) finds that 



where Y is the ratio of the specific heats, is the body 
temperature, T^ is the free-stream temperature and is the 

free-stream Mach number. 

heaving aside the intricacies involved in the choice i 
of a proper Knudsen number, any classification of the rarefied- 
flow regimes is arbitrary to a certain extent. This arbitra- j 
riness is due to the fact that there is no sharp demarcation j 
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between one regime and t be other. An overall classification 
can be made according to the different physical situations 
that may be present in actual flews. Hayes and Probstein 
(1959) have classified the rarefied flows into the following 
regimes : 


(i) 

Fre e-molecule-flow regime 

(ii) 

First-order-collision regime 

(iii) 

Transitional regime 

(iv) 

Fully-merged-layer regime 

(v) 

Incipicnt-merged-layer regime 

(vi) 

Viscous-layer regime 

(vii) 

Vorti city- inter act ion regime 

(viii) 

Boundary-layer regime. 


Unfortunately, experimental studies of the detailed flow 
field of a rarefied gas streaming past a blunt body are not 
available until recently. In 1969, -Ahouse and Bogdonoff have 
proposed the following classification of the rarefied flows 
on the basis of their experimental study of temperature 
profiles: 

(i) Fr 0 e-mo 1 ec ul e ~f 1 ow regime 

(ii) Transitional regime 

(iii) Merged-layer regime 

(iv) Boundary-layer regime. 
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In this new classification, the transitional regime includes 
the first-order-collision regime and transitional regime 
of Hayes and Probstein (1959). The merged-layer regime of 
Ahouse and Bogdonoff encompasses the fully-merged- layer 
regime, the inc ip ient-merged -layer regime and the viscous- 
layer regime of Hayes and Probstein. Differences between the 
new classification of Ahouse and Bogdonoff and the previous 
classification of Hayes and Probstein are primarily regarding 
the limits established for various regimes. Ahouse and Bog- 
donoff find that the onset of merging occurs at a Reynolds 
number approximately five times greater than the estimate 
of Hayes and Probstein. 

In what follows we shall briefly discuss the nature 
of the various flow regimes. 

Por extremely rarefied flows the mean free path is 
much greater than a characteristic body dimension. Under 
a circumstance, the molecules reemitted from the surface do 
not, on the average, collide with the free-streem molecules 
until they are far from the body. Consequently, it is valid 
to neglect the effect of the reemitted particles on the inci- 
dent stream, at least so far as the effects on the body are 
concerned. This is the basic assumption of the free-molecule 
theory. Por this free-molecule flow, the trajectories of the 
molecules and hence the complete flow characteristics can be 
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determined provided that the particle-body interactions ar e 
known. Analyses by the kinetic theory have given a good un- 
derstanding of the aerodynamic characteristics in the free- 
molecule flow. However, calculations depend upon the parti- 
cular model of surface interaction assumed. Expressions for 
heat transfer and aerodynamic forces are obtained in terms 
of overall quantities like the momentum and energy exchange 
accommodation coefficients. Information on the free-molecule 
flow theory may be found in texts such as by Patterson (1956) 
and Schaaf and Ghambre* (1958). 

In the transitional regime between the free-molecule- 
flow regime and the merged-layer regime the mean free path 
is of the same order as a typical body dimension. Surfs'-* ~ 
collisions and free-stream int ermol ecular collisions are of 
equal importance. Air molecules reflected upstream by the body 
now collide with other oncoming molecules, instead of trave- 
lling unhindered. There is a tendency for the molecules to ; 

pile up in front of the body. With increase in the ambient 
density, this pile up increases and the upstream extent of j 

this disturbed region contracts towards the body. The particle 
number density in the pile up region is much greater than j 

in the undisturbed free stream. The portion of the transitional 
regime immediately adjacent to the free-molecule-flow regime, 
lends itself to some special theoretical models. Por example - 
we have the first-order-collision theory of Baker and Charwat 
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(1958), which involves a detailed consideration of collisions 
between certain classes of molecules in the flow field. In 
particular, only the collisions between free-stream molecules 
and molecules directly emitted from the body are considered. 
Molecules reemitted from the surface of the body occasionally 
strike incident molecules, deflecting them either away from 
the surface or into the. surf ace. It happens that the for: 
is the dominant process, producing a shielding effect which 
reduces heat transfer and drag as the density increases above 
frec-molecule limit. Baker and Charwat (1958) derived the 
following analytical expression for the transitional correc- 
tion to the drag coefficient of a sphere at high Mach 
numbers and for small values of the body to free-stream tern- 
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The merged-layer regime is characterised by a 
shock wave and a thick boundary layer. As the free-stream 
density is decreased, the shock wave . increases in thickness 
more rapidly than the boundary layer. Eventually both of them 
merge into each other, and the inviscid region between them 
disappears. As the boundary layer and the shock wave grow In 
thickness, the total extent of the upstream disturbance of the 
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body increases. In this regime the use of neither the Ran- 
kine-Hugoniot relations nor the classicial boundary- layer 
approach is valid. Rarefaction effects like the slip velocity 
and temperature jump conditions are significant in this 
regime. 

It is usual to subdivide the merged-layer regime 
into three parts, namely, the viscous-layer, the incipient- 
merged-layer and the fully-merged- layer regimes. These regimes 
owe their existence more to certain theoretical models that 
have been applied to the problem than to physically demar- 
cated situations. The viscous-layer model is proposed by 
Hayes and Probstein (1959). Here an infinitely thin bow shock 
is supposed to be followed by a fully viscous shock layer. 

The viscous-layer model is the simplest of the models to 
describe the hypersonic viscous flow in the stagnation region 
of a blunt body. It is used by several investigators like 
Probstein and Kemp (i960), Ho and Probstein (1.960) and Shih 
and Krupp (1967)* In the light of the theoretical work by 
Van lyke (1962) and Bush (1964), this model is dropped 
subsequently as being physically inconsistent. In the inci- 
p ient-mer ged-lay ex regime, the shock wave is thick and cannot 
be treated as a discontinuity satisfying the Rankine-Hugoniot 
relations. In one approach to this regime, the shock is 
treated as a discontinuity but the shock relations are modi- ; 
fied taking into account the viscous stresses and the heat i 
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conduction immediately behind the shock wave. In this case 
the total enthalpy and the tangential velocity are not con- 
stant across the shock wave. This regime is studied by 
Probstein and Kemp (i960) and Cheng (1961). In the fully- 
merged -layer regime, the shock thickness is comparable to 
and perhaps even greater than the shock layer thickness* 

Not only is the shock layer fully viscous , but also the shock- 
transition zone may no longer be distinguished from other 
parts of the flow field. 

The boundary-layer regime contains the classical 
boundary-layer regime as well as the vorticity- interaction 
regime. In the classical boundary-layer regime, one finds 
the well-known flow pattern of an infinitesimally thin bow 
shock, a thin boundary layer adjacent to the body and an 
inviscid potential flow in-between. The flow details within 
the shock structure need not be considered and the usual Ran- 
kine-Hugoniot jump conditions for the shock can be used. The 
inviscid and the boundary-layer flows can be computed in a 
piece-meal fashion. A detailed analytical treatment of the 
inviscid hypersonic flow in the stagnation region of a blunt 
body is given by Li and Geiger (1957). Various procedures 
to calculate the inviscid flow are discussed by Hayes and 
Probstein (1959). A thorough and very general investigation 
of the boundary-layer flow for the hypersonic case is made 
by Lees (1956) and Pay and Riddell (1958). In the vorticity- 
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interaction regime, the vorticity generated by the curved 
shock wave in the inviscid part of the shock layer is suffi- 
ciently high and the usual boundary- layer concepts must be 
modified. As a simple model of the phenomenon, Id (1955a, 1955b) 
considers the flow past a semi- infinite flat plate in a 
uniform shear flow, wherein he omits the pressure term in 
the boundary-layer equations# later Id (1956) adds a stream- 
wise pressure gradient which is induced in the boundary layer 
due to the interaction between the displacement thickness 
and external vorticity. However, the existence of this pre- 
sure gradient is disputed by G-lauert (1957). Finally, using 
the method of singular perturbations, Murray(l961 ) and later 
Toomro and Rott (1964) have confirmed the existence of Li's 
induced pressure gradient for an unbounded shear flow. For 
a free-stream shear flow of finite width, the induced pressure 
gradient is negligible. The effect of external vorticity in 
the stagnation region of an axisymmetric Incompressible flow 
is studied by Rott and lenard (1959) and Kemp(l959), They 
find that the presence of external vorticity increases the 
shear at the body surface and that this extra shear is about 
0.649 times the external value of the shear* : 

A diagram given by Ahouse and Bogdonoff (1969) 
depicting the various regimes of rarefied-gas flow, in terms 
of Mach numb er-Reynol ds number axes, is reproduced In 
Figure 1.1, to indicate the general nature of these regimes. 
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This demarcation of the flow regimes is partially guided 
by the cold-wall, temperature-field studies of Ahouse and 
Bogdonoff near the stagnation region of a Blunt Body, in 
the Mach number range of 18-25. Ihile this classification 
suffers from its natural arbitrariness, it need not Be valid 
for other body shapes and other temperatures. Schaaf (1963) 
mentions that a diagram classifying the rarefied-flow regimes 
should atloast Be three dimensional, with some sort of a 
thermal parameter as the third axis, and that even this would 
be entirely inadequate if real-gas effects are to Be included 
The present knowledge of the rarefied gas flows, theo 
Qr experimental, is not exhaustive enough to warrant such a 

classification. 

The investigation of the nerged-layer regime is the 


principal aim of the present work. The various subdivisions 
of the morged-layer regime are icteric! in the present work, 
whore a unified numerical approach, valid for the entire 
mergod-layor regime, will he adopted. The rarefaction effects 

, • +Vii s regime. Here, there are three 

play an important role m this regi 

, , -l,,.!. interrelated parameters, namely, 

separately important But mterreiau 

tho Mach number, the Reynolds number and the appropriate 

Knudsen number. These parameters indicate the importance of 

. . OAaitv an d rarefaction respectively. As 

compressibility, viscosity ana rar 

4- * nol +n 't'VlB 2?3,"fciO of "tllO 

the Knudsen number is proport 10 

number, a moderately laxge Knudsen 

number to the Reynolds number, 



22 


number would mean that either the Mach number is large, 
or the Reynolds number is small or both. Hence the rare- 
faction effects are associated with very strong compressi- 
bility and viscous effects. 

1 Evaluation of the Existing Theories 

Existing continuum approaches to the hypersonic, 
low- lley no Ids -number, blunt-body problem can be divided into 

three main categories. They are 

1 . Higher-order boundary- layer theory 

2 . Two-thin-layer theory 

3. Numerical integrations through the shock. 

A survey of these techniques is made by Cheng (1966) and 
Jain (1968). Some useful discussions can be found in the 
review paper by Mikhalov et al. (1971). In what follows,- a 
critical evaluation of the work done along these three lines 
is given, leading to the motivation of the present investiga- 
tion, 

llikhor-Order Bound ary- layer Theory :- 

The classical boundary- layer theory of Prandtl can 
bo embedded as the leading approximation in a systematic 
asymptotic expansion of the Navier-Stokes equations in an 
inverse power of the Reynolds number . The method of inner 
and outer expansions, also called the matched asymptotic 
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expansions (Van Dyke, 1964a) is used here. Different types 
of expansions are needed in different regions of the flow, 
depending on the order of magnitude of the thickness of these 
regions. If the Reynolds number is not large enough for 
the classical boundary-layer theory to be valid, more than 
one term of the expansion is required. So far, calculations 
arc made upto the third-order boundary layer, treating 
Prandtl *s boundary layer as the first. Exhaustive work on 

the second-order boundary layers is done by Van lyke (1962, 

* 

1963, ',1964b) , Rott and lenard (1962), Lenar d (1962), and 
Mas Ion (1963 ) among several others . Calculations for the 
third-order boundary-layer theory are carried out by Kao 
(1964b) for a cold wall, and Oberai (1964) for an insulated 
wall. Chow and Ting (1961) obtain modified shock conditions 
for piano flow as embedded in the third-order theory. 

Germain and Guiraud (1962) present a more general analysis, 
of the modified shock conditions for three-dimensional flows. 

The procedure followed in the higher-order boun- 
d ary -layer theory may be outlined as follows. The flow 
quantities arc expanded asymptotically for large Reynolds 
numbers in the shock wave, the inviscid region and the 
boundary layer. Appropriately stretched variables are used . 
for the expansions in the shock wave and the boundary layer - 
It is estimated from the previously known first-order solu- 
tions that the thickness of the boundary layer is proportional 
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to 1/lfRe, and the thickness of the shock. wave to 1/Re, .. . 

where Ro is a Reynolds number (von.Mises, 1950). The esqoan- 
sion schemes and the magnifying . coordinates are taken in 
accordance with these estimates. These expansions are sub 
stituted into the Navier-Stokes equations and the like terns 
are equated giving . different sets of equations in different 
regions of the flow. These equations are solved according 

to the following scheme: .. .. 

1. First-order boundary-layer approximation. - a) the 
basic invise id flow is calculated, and b) the. first-order . . 
boundary- layer equations are solved matching with the basic 


inviscid flow. 

2. Second-order boundary-layer approximation:- a) the 
second-order equations in the inviscid region are solved 
taking into account the displacement caused by the thickness 

of the first-order boundary layer, and b) the second orde 

, , match ins with the changed 

boundary-layer equations are solved maten g 

inviscid flow. 

3. Third-order boundary-layer approximation, a) the 

4 -v, o p+micture of the shock wave 

first- order equations for the otr 

are solved, matching with the free stream ahead and the basic 

inviscid flow behind, t) the third-order equations in the , 

4 + into account the modified shock-wave 

inviscid region, taking into aco 

, , -f i -rst— order shock— structure 

conditions which depend on the f 
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solution and the shock curvature, c) the third-order boun- 
dary- layer equations are solved matching with the third- 
order inviscid flow, and finally, d) the third-order equations 
for the structure of the shock wave are solved, 

For any given Reynolds number the above component 
solutions can be combined to give a composite solution that 
is uniformly valid throughout the flow field. 

The higher-order boundary-layer theory can describe 
correctly the flow at moderately small Reynolds numbers, in 
which the interaction between the boundary layer, the invis- 
cid region and the shock wave is small. The theory allows 
identifying several second-order effects and calculating 
separately the influence of these on the flow field. Van 
Dyke (1962) distinguishes seven second-order effects as 
follows: 

1 . Longitudinal curvature, 

2. Transverse curvature, 

'5 • Di,3 p 1 ac cm ent , 

4, Ext crnal gradient of entropy, 

5 * External gradient of enthalpy , 

Slip velocity, 

7 . T omp erature j ump . 

External vert lefty is related to the external gradients 
of entropy and enthalpy through the Crocco ' s relation. These 
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effects (except the slip velocity and temperature jump) 
are not independent from each other in a strict sense. 

The difficulties resulting from the arbitrary division of 
these effects can be avoided by evaluating then concurrently. 
Van Dyke (1963) gives the following expression for the sta- 
gnation-region heat-transfer rate for a sphere. at = 4 
in a perfect gas withy= 7/5, Prandtl number 0.7, viscosity 
proportional to temperature and the surface cooled to 0.5 
times the f roe-stream stagnation temperature: 


external curvature slip and displ. 
vorticity temp. jump speed 


— SL. 

q BL 


= 1 + [ 0.57 


+ 0.17 - 0.39 - 0.18 


- 1/2 

0.17]Re g 


If the body is cooled instead to 0.2 times the stagnation 
temperature, this expression is replaced by 

— 1/2 

— » 1 + [ 0.46 + 0.08 - 0.23 + 0.00 = 0.31 ] Re , 

q BL 

where do., is the Reynolds number based on conditions behind 
the: normal shock and the body radius. External vorticity 
is soon to be the leading effect. A detailed discussion of 
the various second-order effects can be found in the review 
papers by Chong (1966), Jain (1968) and Van Ryke (1969). 
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The higher-order boundary-layer theory does not 
give accurate results for low-Reynolds-number flows corres- 
ponding to the fully-nerged-layer regime. The essential 
feature of this regime is that the entire flow field is 
viscous. The boundary- layer theories envisage an essentially 
inviscid zone between the shock wave and the boundary layer, 
oven though in the third-order boundary- layer theory viscous 
terras appear in this zone. In the fully-merged-layer regime 
the departure of the flow quantities from the classical 
boundary- layer solutions are large and cannot be treated as 
perturbations. Further, the method of treating various zones 
of the flow piecewise is not valid in this regime. Indeed, 

Kao (1964b) finds that for Re g less than about 100 the sur- 
face heat-transfer rates predicted by the higher- order 
boundary-layer theory deviates remarkably from the experi- 
mental values. Recently, Davis (1972) numerically integrates 
the full Navicr^Stokes equations for the incompressible flow 
past a parabola. He compares his results with Van Dyke *s 
(1964b) second-order boundary layer results and finds that 
the later diverge even at Re^ - 1000. The failure of the 
higher-order boundary- layer theory to describe the low Rey- 
nolds number flow expresses itself in an . interesting way. 

Such anomalous situations occur as the boundary-layer over- 
taking the shock position or the thick shock wave encompass- 
ing even the body (e.g. Kao, 1964b). In such cases meaningful 
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composite solutions cannot be obtained. 

Two-Thin- Layer Theory 

Notable investigations using this model are carried 
out by Cheng' (1961, 1963a, 1966), Chow (1963) and Liu (1967). 
All these analyses use the continuum equations of Navier- 
Stokes with a thin-layer assumption, which greatly simplifies 
the governing equations. The main feature of this theory 
is to divide the flow into two adjoining regions, namely, 
the shock wave (or the shock-transition zone) and the shock 
layer, both of which are assumed to be thin. In the shock- 
transition zone it is assumed that the mass, momentum and 
energy fluxes tangential to the shock wave are neglected. 

This assumption, together with the thin-layer approximation, 
reduces the governing equations in the shock-transition zone 
to a set of ordinary differential equations similar to that 
for a plane shock. The governing equations for the shock 
layer, that is, the region between the inner edge of the 
shock wave and the body, embody differential equations of 
both the inviscid shock-layer theory and the classical 
boundary-layer theory. These are obtained using the thin- 
layer approximation and a high, shock-compression ratio. These 
equations still remain parabolic and their mathematical 
analysis is reduced to a problem compared to that encountered 
in the classical boundary-layer theory. 
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An essential feature that the two-thin-layer 
theory incorporates, and the higher-order boundary-layer 
theory fails to do properly, is to preserve' the viscous 
and energy transport effects just behind the shock. The 
shock-transition-zone equations lead to a modified form 
of the Rankine-Hugoniot shock conditions as follows: 

P 8 V S = V co 

p s = p oo vi 

^OO V 00 = ( U JjJyOg 

p <» (H s - HJ =( §*: 5 y [ H - (1 - Pr) 4 ] > s 

where the subscripts 00 and s refer to the conditions at 
the free stream and just behind the shock-transition zone 
respectively, y is the length normal to the shock, and H 
the total enthalpy. According to these modified relations, 
the tangential velocity and the total enthalpy are not con- 
served across the shock but are controlled by the transport 
effects behind the shock which are not generally small. These 
terms being proportional to the velocity and temperature 
gradients as in the case of slip boundary conditions, the 
modified Rankine-Hugoniot relations are sometimes interpreted 
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as shock-slip conditions (Pan and Probstein, 1962). 

It is interesting to note a simplification affor- 
ded by the two -thin-layer theory. The flow in the shock 
layer together with the location of its outer edge can be 
determined entirely independent of the flow in the shock- 
transition zone. The modified Rankin e-Hugoniot conditions 
and the appropriate conditions at the body surface (for 
example, the vanishing of the normal velocity, and the slip 
boundary conditions) are enough to solve the shock-layer 
equations completely. Nevertheless, the flow in the shock- 
transition zone depends on the shock-layer solution. The 
later provides explicit boundary conditions to integrate the 
Ordinary differential equations that govern the shock-tran- 
sition zone. 

With this basic formulation of the two-thin-layer 
theory Cheng (1961) investigates analytically the shock- 
layer flow in the stagnation region in the incipi ent-merged- I 
layer regime. It is found that the slip-like boundary condition! 
at the shock, when used together with no-slip boundary condi- 
tions at the body surface, lead to a unique property of the I 

solution, namely, in the limit of zero Reynolds number, the ! 

surface heat-transfer and skin-friction characteristics 
tend to the correct free molecule values for unit accommodatioi 
and reflection coefficients. Further Cheng and Chang (1964) { 

give an analytic solution for the shock-transition zone in 
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the Newtonian (Y -* 1) and the hypersonic (M^ -* °°) limits* 

This analytic solution is compared with Levinsky and Yoshi- 
hara’s (1962) direct-integrat ion- thro ugh- the-shock results, 
which are obtained for a monatomic gas (y =1,667). limit- 
ed difference between the two solutions is somewhat surprising 
considering that y = 1.667 is rather far from unity, near 
which Cheng *s analysis is supposed to be valid. Cheng and 
Chang (1964) also study the effect of slip boundary conditions 
as correction terms in a perturbation scheme. They obtain 
the following expression for the slip correction to the heat- 
transfer rate i 
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where Ctt is the no-slip value of the heat-transfer coeffi- 
n o 

cient given by 
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X and are the proportionality constants in the expre- 
ssions for slip velocity and temperature jurnp^ and the j 
mean free path corresponding to the temperature and pressure^., 
at the body surface. |l/3 l (k) is the incomplete gamma function 
of order 1/3 in the argument k. (Other symbols are defined 
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in the nomenclature.) liu (1967) investigates the effect 
of wall temperature on the flow in the shock layer within 
the framework of Cheng’s theory. Davis (1970) studies the 
shock layer essentially with Cheng’s formulation, hut with 
more general governing equations than used by Cheng, Both 
Diu and Davis incorporate the slip boundary condition effects 
as higher-order terms in a perturbation scheme. 

Bush (1964) formulates the hypersonic viscous blunt - 
body problem with asymptotic expansions of the variables in 
terms of 1/Re, 1 /M^ and To make the strict asymptotic 

expansions possible, Bush divides the shock-transition zone 
into three subregions, namely, (A) the outer region, (B) the 
middle region and (c) the inner region. Including the shock- 
layer region, Bush thus analyses the flow on the basis of 
four distinct regions as opposed to only two in Cheng’s for- 
mulation of the two -thin-layer theory, Cheng (1966) finds 
that his equations of the shock transition zone amount to 
’compo.site equations', which are uniformly valid in the 
regions (A), (B) and (C) of Bush. Thus, Cheng establishes the 
equivalence of his two-thin-layer theory with the formulation 
of Bush, 

The two-thin-layer theory, with its simple . formulation 
and amenability to analytical form of solution in several 
cases, predicts the stagnation point heat -transfer rates 
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in good agreement with experiments, and is in several 
respects superior to the higher order "boundary- layer theory. 
However, one should not expect the two-thin-layer theory 
to predict the detailed flow field correctly at least in 
the fully-merged -Hay er regime, because of several simpli- 
fying factors which are not valid under rarefied conditions 
of flow. In particular, we mention the following: 

(1 ) The basic assumption of a thin layer of distur- 
bance is not valid under very rarefied conditions, where the 
shock wave and shock layer are very thick and are comparable 
to the body nose radius* 

(2) The shock-thickness effects and the resulting 
effects of tangential fluxes of mass, momentum, and energy 
across the shock-transition zone are neglected. 

(3) The coupling between the solutions in the shock- 
transition zone and the shock layer is weak. The flow in the 
shock layer can be solved independent of the shock-transition 
zone, After calculating the flow in the shock layer, the stru- 
cture of the shock transition zone is studied. In most of 

the calculations there lies a discontinuity in temperature 
and normal velocity at the interface of the shock-transition 
zone and the shock layer. When the flow is fully merged, one 
expects that the shock structure should exercise significant ; 
effect or the shock layer. 
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(4) All the investigations using this model of flow 
take slip velocity and tempqn^turo jump as small. and evalu- 
ate corrections in the flow due to these effects, This may 
be justified only for a cold body> 

The failure of Cheng's two-thin-layer model to 
predict the correct behaviour of the entire flow field is 
exemplified by a comparison with recent experiments reported 
by Ahouse and Bogdonoff (1969). In Figure 1,2, reproduced 
from Ahouse and Bogdonoff *s paper, it can be seen that 
Cheng's theory predicts far too thinner a profile compared 
to the experiment, Ahouse and Bogdonoff have attributed 
this to the failure of the Navier-Stokes approach to describe 
the morged-layer flows. But this conclusion may not be 
correct because Cheng's theory does not represent the full 
Navier-Stokes equations. 

Numerical Integration Through the Shock :- 

In this approach no special assumptions are 
needed regarding the model of flow. The governing equations, 
either the full Navier-Stokes equations or some simplifi- 
cations thereof, are integrated numerically from the body 
to the freest ream, where uniform conditions exist. Depending 
on the flow parameters like the Mach number and Reynolds 
number, the position of a sharp shock or the merged profile 



of a thick shock wave and thick viscous layer near the 
body come as an integral part of the numerical solution 
itself. The difficulties associated with the numerical 
integration of the full Navier-Stokes equations are prohi- 
bitive. As such, most of the work using this approach is 
limited to a region near the stagnation point. Here the. 
assumption of local similarity reduces the Navier-Stokes . 
equations to a system of ordinary differential equations. 

The investigation of only the stagnation region does not 
impose any highly serious restriction to the basic understa- 
nding of the rarefaction effects in the flow field. The 
stagnation region, where the heat transfer is maximum, is 
the most important region of a blunt-nosed reentry vehicle. 
Also, the stagnation region studies provide initial solutions 
for integrating the governing equations donwstream. Moreover, 
most of the important experimental results for the rarefied, 
blunt-body problem are available only for the stagnation 
region. 

So far the approach of numerical integration 
through the shock is the most reliable approach to a theore- 
tical understanding of merging and related rarefaction eff- 
ects. However, the price of this reliability is paid by the 
complexity in the numerical integration and by the large 
amount of computer time needed, levinsky and Yoshihara (1962) 
and Kao (l 964b) use this approach in their investigation of 



36 


the blunt -body problem. Also Chung et al. (1968) and Delli- 
nger (1971) use this approach to study the non-equilibrium, 
merged-layer, stagnation-region flows. While the numerical 
techniques used by these authors differ, the unifying factor 
is that the entire flow field is determined by integrating 
the same system of differential equations. 

Lev insky and Yoshihara (1962) obtain exact numer- 
ical solutions of a system of simplified Ravier-Stokes equa- 
tions for insulated and cold wall cases in a monatomic gas 
e ream at Re^ = 152, 1382, 13682 and = 10. These are the 
first solutions obtained using the method of numerical inte- 
gration through the shock and have for a long time served 
as a basis for assessing the adequacy of approximate theo- 
retical models (sec, for example, Oberai 1964, Cheng 1966). 
Levinsky and Yoshihara use a thin-layer assumption to reduce 
the Ravier-Stokes equations to a simpler form. Nevertheless, 
the thin-layer equations of Levinsky and Yoshihara include 
many terms of the original equations that Cheng* s two-thin- 
layer model neglects. This extends the validity of Levinsky 
and Yoshihara* s theory to much lower Reynolds numbers than 
Cheng *s theory in predicting the detailed flow structure. 

At Re TO s 152 one finds that, for the insulated-wall case, 
the region of disturbance is as much as about 50 percent of 
the body radius. Further integrations of Levinsky and Yoshi- 
hara *s equations at lower Reynolds numbers, made as a 
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complementary part of the present investigation, show that 
there is a steep increase of the extent of upstream distur- 
bance. Thus the thin-layer assumption used by Levinsky and 
Yoshihara is of doubtful validity when Re^ is less than 150. 

A further confirmation of the inability of Levi ns ky and 
Yoshihara* s theory to predict the very low Reynolds number 
flow cones from a recent investigation of Vogenitz and Takata 
(1970). They calculate the merged-layer flow-field with Monte 
Carlo simulation technique, allowing for intermolecular col- 
lisions. They find that the Monte Carlo profiles differ very- 
much from those of Levinsky and Yoshihara, and attribute 
this difference to the failure of the continuum model to 
describe the merged-layer flow. As in the case of Ahouse and 
Bogdonoff *s (1969) temperature-profile comparisons, this 
conclusion nay not be true because Levinsky and Yoshihara 
do not use the full Navier-Stokes equations. Further, it 
may be noted that Levinsky and Yoshihara's results predict 
a consistent drop of the wall pressure as the Reynolds number 
is decreased, both in the insulated and the cold wall cases. 
This is contrary to the experimental trend established by 
Sherman (1953), and later by Potter and Bailey (1965) and 
several others. Levinsky and Yoshihara have conjectured that 
this discrepancy with experiments, may be due to the neglect 
of slip effects. As a result of the present investigation 
we know that this is due to the thin- layer approximation 
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rather than the neglect of slip velocity and temperature 

jump, 

Kao (1964b) uses the full Navier-Stokes equations 
together with the usual local-similarity assumption in his 
numerical integrations through the shock. He obtains. solutions 
only for three representative Reynolds numbers,. Re = 10,100 
rand 1000, where Re g is based on conditions behind a Rankine- 
Hugoniot shock and nose radius, at = 10 and y= 11/9. 
Moreover, Kao considers only the case of a very cold wall, 
and does not include slip boundary conditions. These compu- 
tations are not exhaustive enough to give an understanding 
of the merged flows. 

We observe that the work so far carried out using 
the method of complete numerical integration, suffers from 

the following deficiencies : 

(1) The numerical solutions are obtained only for few, 
stray values of the parameters like the Reynolds number and 
the Mach number. They are not exhaustive enough to understand 
fully the mechanism of the raerged-layer flows. 

(2) The solutions neglect the slip velocity and tem- 
perature jump. This is an unfortunate omission, particularly 
considering the potentiality of this method to allow these 
boundary conditions as principal terms. 
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(3) The effect of the wall-temperature variation on 
the flow is not studied, 

(4) Solutions obtained using the thin-layer assumption 
are not valid at very rarefied flow conditions, where the 
upstream extent of disturbance is comparable to the body 
nose radius. 

( 5 ) Solutions obtained using the thin-layer assumption 
predict a drop in the surface pressure as the Reynolds 
number is decreased. Experiments show first a decrease and 
then a steep increase over the high Reynolds number limit. 

1 .5 Survey of the Experimental Work 

Experimental results for the hypersonic, rarefied, 
blunt -body problem are available for the heat-transfer rate 
to the body, the surface pressure and the detailed flow- 
field studies like density and temperature profiles. Prom 
literature survey, we find that no detailed flow-field stu- 
dies are made until 1967. With the development of electron- 
beam, X-ray and fluorescence techniques, .an entirely new 
approach to problem has become available, allowing detailed 
probing of the flow field. Results for the heat-transfer 
rate and pressure are reviewed by Cheng (1963b) and Potter 
( 1967 ). 
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Heat -transfer rate :- Several experimentalists have 
measured the heat-transfer rate to hemispherical bodies 
in the stagnation region, for various degrees of rarefaction. 
They include Ferri, Zakkay and Ting (1961, 1962), Hickman 
and Giedt (1962), WittLiff and Wilson (1962), Yidal and 
Wit&iff (1963), Valensi and Rebont (1963), Potter and 
Miller (1963), Carden (1966) and more recently Boylan (1971). 
To summarise the situation, we reproduce here , in Figure 1.3, 
a diagram given in the review paper by Potter (1967). Regar- 
dless of the scatter of the data, it may be said that the 
essential features of the aerodynamic heating in the merged- 
layer regime are reasonably clear. As the Reynolds number 
decreases, first increases slightly from the boundary- 

3-BL 

layer value of unity and then decreases towards the free- 
molecule limit. 

Pressure :- Extensive experimental data on the pressures 
in the stagnation regions of hemisphereical and flat-nosed 
axisymetric bodies in rarefied hypersonic flow are presen- 
ted by Potter end Bailey (1963). They also compare their 
results with earlier experimental data. Diatomic and mona- 
tomic gases are used in the study and the experiments are 
conducted with both cold wall and adiabatic wall conditions, 
in the free-stream Mach number range of 4 to 7 and for 
various Reynolds numbers.- 
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Figure 1.4 (from Potter and Bailey, 1963) shows 

the general behaviour of the impact pressure on hemispherical 

noses in rarefied, hypersonic flows. The measured impact 

pressure p^ is nondimensionalised with respect to the invis- 

cid impact pressure p* and is correlated in terms of the 

, p sJ/ § 

■.parameter Re (—r— ) , where p o is the density behind the ■ 

° M oo b 

shock wave and Re g is the Reynolds number based on the nose 
radius and conditions immediately downstream of a Rankine- . 
Hugoniot shock* The most important feature of the data is 
that, as the Reynolds number is decreased, the measured impact 
pressure first decreases slightly from the inviscid limit 
and then begins to rise rather rapidly. This rise occurs at 

, Pc 1/2 

Re s (— -~) mi 100. Potter (1967) suggests that the initial 
decrease in pressure is due to shock weakening. Approximate 
analyses like those of Probstein and Kemp (i960), levinsky 
and Yoshihara (1962) and Oberai (1964), which take into 
account the transport effects behind the shock, lead to a 
decrease of the impact pressure in the incipient-merged- layer j 
regime. However, they fail to predict an increase of the j 

impact pressure as Re g is further decreased and the fully- j 

merged-layer regime is approached. Thus one concludes that j 

in the fully-merged-layer regime, rarefaction phenomenon f 

affects the mechanics of the flow rather intrinsically, whichj 
is not predictable by the simplified analyses. Probably, in ; 
this regime the structure of the shock plays an important j 
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role in determining pressure. In the present investigation, 
we find that the pressure first decreases and then increases 
rapidly as the Reynolds number is decreased. 

Detailed flow-field measurements Electron-beam probes 
are used by Oguchi et al. (1965) and Ivanov (1965) to study 
the density field in the shock layer of blunt bodies, each 
of them for a single value of Reynolds number and Mach number. 
These investigators measure the absorption of the beam. It 
involves integrations of the data along the beam path to 
yield the local density. A more complete and comprehensive 
study of the density field is reported by Wainwright (1967), 
who measure the fluorescence created by the beam. The fluo- 
rescence technique does give a local measurement of the 
density, but the cloud of secondary electrons surrounding 
the main beam can cause an adverse influence, 

Russel (1968) uses an 2L-ray technique instead of 
a fluorescence technique and eliminates the effect of secon- 
dary electrons, Russel gives the detailed merged-layer density! 
profiles for argon and nitrogen along the stagnation stream 
line of spherical models. He uses the nitrogen stream with 
a nominal Mach number of 4.2 and the argon stream with 3.8. 
Density profiles are provided for a Re^ range of 30-500 in 
nitrogen and 100-2000 in argon. The measurements are made 
with model surface temperatures of 300°K (equal to the free- 
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stream stagnation temperature), and of 78°K. In Chapters 4 
and 5 we compare our profiles in detail with those of Russel. 

Ahouse and Bogdonoff (1969.) measure the local 
rotational and adiabatic recovery temperatures along the 
stagnation stream line of a hemispherical nose in a nitrogen 
stream, using the electron beam fluorescence technique and 
a radiating hot wire probe. In most of the measurements a 
free-stream stagnation temperature of 1600°K is maintained. 
Variation of the free-stream stagnation pressure gives a 
Re 0 variation of 50 to 460 and a BC variation of 18 to 26. 

The wall temperature is kept at about 0,19 times the free- 
stream stagnation temperature. Detailed comparison of our 
results with these profiles is presented in Chapter 4. 

Attempts to extend Ahouse and Bogdonoff *s work to 
still lower Reynolds numbers are recently reported by Becker 
et al. (1970). from this review of literature, we. find that 
there is still much dearth of experimental data to give reli- 
able information about the detailed flow structure at low 
Reynolds numbers. 

1.6 Statement of the Brohlem and Scope of the Present Work 

The flow in the stagnation region of a sphere, 
placed in a rarefied hypersonic stream. of a perfect gas, is 
investigated. In the previous sections, several limitations 
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of the existing continuum studies of the rarefied-blunt - 
body problem are pointed out. Some of the significant 
questions that still await answers are 

(1) What is the limit of validity of the Navier- 
Stokes equations in the hypersonic, rarefied flow around 
a blunt body? 

(2) How far can the thin-layer assumption, that is 

so widely used in the raref ied-blunt-bcdy problem, be expec- 
ted to give reasonably good results? 

and (3) What is the effect of the slip velocity and 
temperature jump boundary conditions on the flow field? 

The existing theories are unable to give satisfac- 
tory answers to these questions owing to two reasons. 

Firstly, they postulate several simplifying models which do 
not correctly represent the physical situations existing in 
the merged-layer flow. Secondly, the slip boundary conditions 
are either neglected or assumed to be small and corrections 
are sought to the no-slip solutions, instead of treating 
them as principal boundary conditions. 

The present investigation is motivated by the 
search to find answers to the questions posed above. The 
full Navier-Stokes equations are assumed to describe the 
hypersonic , rarefied flow past a sphere. The region of 
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interest is the stagnation region where the properly non- 
dimensional is ed governing equations can he reduced, using 
local-similarity concept (Kao, 1964a), to an eighth-order 
system of non-linear coupled ordinary differential equations, 
with two-point boundary conditions. Slip boundary conditions 
are used at the surface of the body. Either the body temper- 
ature is prescribed or it is given that the body is adiaba- 
tic. At the other end, the free-stream infinity, the flow 
quantities asymptotically obtain their uniform free-stream 
values. Two important parameters that are to be prescribed 
are the Mach number and the Reynolds number. Numerical solu- 
tions are sought neglecting the slip boundary conditions as 
well, to know the effect of slip on the flow field. These 
solutions will be usually referred to as the no-slip solutions* 
In addition, simplified equations are derived using a thin- 
layer assumption in the manner of Levinsky and Yoshihara 
(1962), smd complementary numerical integrations are carried 
out for these equations with and without slip boundary 
conditions. Comparison of the solutions of the full Navier- 
Stokes equations with solution of the thin-layer equations 

will determine the limit of validity of the later. | 

I 

In our investigations we consider only a perfect 
gas with constant specific heats. A more realistic analysis 
of the low-density hypersonic flow must include nonequili- j 
brium flow chemistry with various dissociation, ionisation 
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and recombination processes of the air molecules. However, 
the effect of these nonequilibrium processes on the flow are 
only quantitative in nature and the fluid-dynamic aspects 
are not intrinsically changed. The aim of the present inves- 
tigation is a basic understanding of the mechanics of the 
merged-layer flows, and we restrict ourselves, as a first 
approximation, to a perfect gas with constant specific heats. 
Some investigations of the real gas . effects on the blunt- 
body flows are made by Chung (1961), Cheng (1963),. Shih 
and Krupp (1967), Chung et al. (1968) and Dellinger (1971 )• 

As a result of the present investigation we find 
that the slip boundary conditions affect the entire structure 
of the flow field (including the shock-wave-like region). 

These effects increase quantitatively as the Reynolds number j 
is decreased. In the adiabatic -wall case, slip velocity j 

increases to as much as 60 percent of the free-stream 
velocity. Slip velocity and temperature jump have opposing 
effects on the flow field. For example, the slip velocity j 
has a tendency to lessen the total extent of disturbance, 
while the temperature jump increases it. .The total effect j 
of the slip boundary conditions depend upon the wall tempe- j 
rature. Compared to a cold-wall case, slip effects are more 
on an adiabatic wall. 

It is found that the thin-layer and the two-layer-., j 
model theories fail to describe correctly the detailed flow j 
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field at low Reynolds numbers, even though they are able 
to predict certain surface characteristics like the heat- 
transfer rate to the body in good agreement with experiments. 
Detailed flow comparisons .with the present exact solutions 
show that the solutions using the thin-layer assumption are 
not valid below Re TO = 150, Higher-order boundary- layer 
theories are found inadequate to describe the flow in such 
a low -Reynolds -number regime. 

Comparison of the present results with merged-layer 
experimental density profiles of Russel (1968) and tempera- 
ture profiles of Ahouse and Bogdonoff (1969) show a good 
agreement, The present results also show a steep rise in the 
impact pressures and the adiabatic recovery temperatures at 
low Reynolds numbers as observed in experiments (e.g. Potter 

and Bailey, 1963, and Drake and Backer, 1952). The thin-layer ; 

'■ | 

theories fail to predict such features. even qualitatively. 

These factors establish the superiority of the full Navier- 

l 

Stokes equations over the thin- layer .theories in describing 

' "■ ' ' ■- ' ■ ■ \ 

the mechanics of a fully-merged flow. 

In Chapter 2, we present the mathematical formulation 
of the problem. In Chapter 3, we discuss the numerical 
technique employed to integrate the differential equations. 
Since we aim to obtain extensive solutions for wide ranges ! 
of Reynolds number, Mach number and other parameters, the 



proper choice of an efficient numerical technique is of 
paramount importance, in view of the computer time needed* 

The main numerical difficulties arise from the fact that, 

1 ) 'the equations are nonlinear and coupled, 2) the boundary 
.conditions are given at two points, one of which is an 
unknown 'free-stream infinity*, where all variables should 
tend asymptotically to their free-stream values, 3) the 
slip boundary conditions are given in a differential form 
and are a priori unknown, and 4) steep gradients as in the 
case of a thin shock wave give rise to numerical instabili- 
ties. We employ a finite-difference relaxation scheme, re- 
cently coming into prominence (e.g. Dellinger, 1971). This 
method, known as the Accelerated Successive Replacement 
method, is most perfectly suited for boundary-value problems 
of present complexity. It has the following salient features: 

1 ) in each individual iteration, the coupling of the equations 
is eliminated, 2) the boundary conditions are always satis- 
fied, 3) the corrections applied in each iteration are con- 
trolled to guard against any blow-up owing to numerical 

i 

instability, and 4) it is easy to obtain the asymptotic 
nature at the free-stream. The simplicity of this method 
and the gain in the computer time are amazing compared to . 
other methods that can be possibly applied* In Chapter 3 f 
in addition to describing the numerical technique and keeping 
the governing equations in a form suitable for the applicatio 
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of this technique, we also discuss some pertinent points 
about the convergence of the scheme. 

The results of the numerical integrations obtained 
for cold-wall and adiabatic-wall cases are presented in 
Chapters 4 and 5 respectively. These results encompass a 
wide range of the Reynolds number, reaching the boundary- 
layer regime on one side and the transitional regime on the 
other. Detailed flow profiles as well as body characteristics 
such as skin friction, heat-transfer rate and surface pres- 
sure are presented. The effect of the slip boundary condi- 
tions and the validity of the thin-layer equations are dis- 
cussed. Extensive comparisons are made with the results of 
other theories and experiments* 

In Chapter 6, the important conclusions drawn from 
the present investigation are reiterated and summed up. 

Thus, this chapter focusses the contribution of the present 
study to the hypersonic, rarefied, blunt-body problem. A 
brief statement on the further research that can be carried 
out is also included. 
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CHAPTER 2 

MATHEMATICAL EORMULATIOI OE THE PROBLEM 


2 , 1 Basic Governing Equations of Motion 

Equations governing the motion of a perfect gas on 
an aocisymmetric "body are the following; 

Continuity equation: 

( p v) r + + ” ( p u) 0 + P u = 0 , (2,1) 

Radial momentum equation: 

2 


= 2 [ yv p - ^ (rv r + 2v + u Q + u cote) 3^ 

1 , r /Uv X * LL. (Orr - Uq + ^ > 

"V* 

+ ~ {cote [r (~) r + -|] - | (v + u cote) } , 

( 2 , 2 ) 


Transverse momentum equation: 



p ( v u r + r U e + 


UV x 

r ■ 



= I 1 1 Cu e + v) - ?r 


i_ ( r , + 2v + u e + u 


cote)] 


/U\ . JL 


lv[r(=) r + ~r 5 


+ .2t tc©, + r-1 


+ cote (u e - u cote), 

T 


(2.3) 


Energy equa 


t ion: 


p (vh r + ^H e ) 


+ r__ u»t n r L+ 


vp r + - Pq + r Pr 


, cote ] + „{ 2 Tj + ^ ( u e + 

+ vlL. + vH e -”1T J + 4 r 

* 2 

v A 


+ -4 (v + u 


cote ) 2 + t r 1 } 


2 v ( r 2 v 2 + 2v 2 + u 0 + u cot 


2 + u 2 cot 2 e + 4 -r v r -v 


+ 2r v r u 


cote + 2 r v r u 0 + u e 


+ 4ut 


cote + 2 u e u cote ^ 


(2.4) 


and the perf 


e ct gas equation 


p = G p 1-7 


(" L ) 


(2.5) 


. 9 and the Telocity components 

The polar coordinates r and , a „ aenote 


• ,^ jre 2 . 1 . Suffices r aim * — 

U and V are shown m +n r and 6 respectively 


... respect to r ana o 

partial differentiation with r ^ ^ 4efto ed 

All the other symbols appear g 
in the list of Symbols. 
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We ma }ce the following assumptions: 


(1) u = v C T )’ 


(2) 0 p 5 y and Pr are constants 


( 2 . 6 ) 


(■ 5 ) Bulk viscosity is zero. 


2 •ppiim^ar-y Oonditinns, 


( 1 ) 


from the body a uniform free stream, with a 


veloc 


y^ towards the body, is a^> 


surned . Thus , as r 


■n = Y sin9 

CO 


y eos0 


(2.7) 


p = p ro , P - p c 


1 = lr 


+ +-hp surface of the hody slip 
(2) (a) 5or the gas at t ^ e rl 
relocity and temperature jump boundary condxtrons are used. 


Thus , at r - m-g 


2 - q x 9u , 

— “ A w 9r 


( 2 . 8 ) 


2 - a 


: w 3T 


TTT Pr 3r 


(2.9) 


,„ here0 end a are reflection and thermal accomodation eoef: 
' c ; 0 nt 3 , and , w is the mean free path of the g as near 

body surface given by 
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A 

W 


V 


TCY 

— 77- 


p~a • 


(b) Further, either the body is given to be adiabatic 
or its temperature, 3^, is prescribed. In the adiabatic wall 
case, 


ar 

9r 


0 . 


( 2 . 10 ) 


Here, the temperature jump in equation (2.9) vanishes and 
the body as well as the gas adjacent to the body have the 
same temperature. In case the wall temperature is prescribed, 

I = prescribed value (2.11) 

w 

(c) lastly, we have the condition that the normal 
velocity vanishes at the wall. Thus, at r = r-g, 

v = 0. (2.12) 

2 . 3 Similarity Postulates and CTondimensionalisation 

For the solution near the axis of symmetry the 
variables are expanded in a series about the axcis of symmetry 
and only terras corresponding to truncation as shown below are 
finally preserved in the equations.. 
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u = u Q (r) sin 8 ' + u 2 (r) 

! 

v = v 0 (r ) cos0 1 + v 2 (r) 

! 

T = T Q (r) J + ^2 ( r ) 

h = h Q (r) * v h 2 (r) 

! 

u = v 0 (r) ’ + u 2 (r) 

f 

p = P Q (r) j + P 2 (r) 


P = P 0 (r) + p 2 (r)sir^' + P 4 (r) 

t 

t 

t 

Truncation 


3 

sm 0 + . . , 

. p 

cos© sm 0 + • . • 

• 2 a 

sm 0 + . . . 

• 2 d 

sm 0 + ... 

- 2 q 

sm 0 + ... 

2 

sin 0 + ... 

• 4 0 

sm 0 + . . . 


( 2 . 15 ) 


Here all the coefficients with subscripts are functions 
of r only. The truncation is equivalent to assuming that the 
flow is locally similar in the stagnation region. 

We make the variables nondimensional as follows: 


v 


u o - Po 


v = 7-5 U= r , P 

v oo 7 CO 


p oo 


T~ ’ 


p,, 

(L 


■o ~ ■? » 

d p Y 

oo v oo 


_ P T 

P = — , I = -j 2 - , 

P oo °°° 


( 2 . 14 ) 


Y x P 

u 0 - r tv- co p oo 


1 v ’ r ’ 

o°° B 


o°° 
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Here, is the adiabatic stagnation temperature 

corresponding to the free stream, and should not be confused 
with.f 0 C°°) which is the value of the function T Q in Equation 
(2.13) as r is the viscosity measured at the tempe- 

rature I . Obviously, a bar denotes a nondimensional variable 


2.4 Full Havier-Stokes Equations 

Introducing the similarity postulates (2.13) and 
the nondimensional is at ion (2.14) into the basic governing 
equations (2.1) to (2.6), and equating the coefficients of 
like powers of 9, we obtain the following system of ordinary 
differential equations, where a prime denotes differentiation 
with respect to r. 


r ( p v) + 2 p (u + v) = 0. 


-r 


U 


— .. 3 Re /— 1 — ~ ' \ /2 u \ - ’ 

v " = - — (p + p v v ) - (- + mr) ^ ■ 

4-J 1 - y 2r 


t 


—i 


+ ( 1 + JL.) u + v 

2r u r 


— f — * p U — \ 

P + p = -tr- (u + vj 
2 r 


u 


u 



u + v 



+ 



V 


u 4- V + 

MW 

r 


ReP u (u + v) 


— •— 

u r 


(2.15) 


( 2 . 16 ) 

( 2 . 17 ) 
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+ B g . p_j . ' v + _£ _ s - c | + C ) + fhl! , 

V 3r r v x y. 


(2.18) 


A P v r 


-_« 2 -2 - _ : 

2 _ ■ * - 2 . 2 v t (v+ u) 

t? = v p r + ~ — + *- 


Re 


Re 


- [ r v + 2(u+v) ] 


+ 


A 


_ t 


Re Rr 


_ 1 - 1 ~ - ~ ti ) 

•V* , , ff’ / J 


r (2vT + r y 1 + r p ! 


(2.1S) 


p = -if-!- iil. 


(2.20) 


Viscosity -temperature relation 


v = V (?) (2.21) 

is prescribed. 

Here 

C T 

A = (2.22) 

Voo 

’A’ can be ea.sily expressed in terms of the free stream 
Mach number as 


ii ^ 



1 

Or-1) iS 


(2.23) 
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The set of equations (2.15) to (2.21) are the main 
equations of the present investigation, and will henceforth 
be referred to as the ’full Navier-Stokes equations 1 . These 
are seven equations to determine the seven unknown variables 
u> v, T,v, p , p and p£* Equations (2.20) and (2.21 ) are al- 
gebraic and they determine p and y . The remaining five 
equations form a set of eighth order, nonlinear, coupled 
ordinary differential equations. 

2 . 5 Boundary Conditions in Nondimensional Form 

In terms of the nondimensional variables the 
boundary conditions (2,7) to (2,12) can be expressed as 
follows : 


(1) 

At the free stream infinity, that is, 

as r ** “ , 

U 5= 

1, 

- (2.24) 

V =' 


(2.25) 

P 2 = 

= o, 

(2.26) 

p = 

= 1, 

(2.27) 

5 = 

1 - 2l“ 1 

(2.28) 

*^t 

y 

1 

o * 

(2.29) 




V = 

v (T„) . 

(2.50) 



63 


(2) On the surface of the body, that is, at r = 1» 
slip velocity; 


5(1) 


2-0 ^ 1 ^ , y d| ) ^ C2v31) 

°°° p - vt dr ^ =1 


Ei/r 


t omp eratur e 3 ump t 


1(1) - 5 W 


r. r 2-a 2 Y 

— 


ML 


/I. 


a Y+1 Pr Re v 1 


QO° 


c-i- i-i 

rjs d? 1=1 

(2.32) 


where, 

r 


w 


T 

x w 

0 00 


For an adiabatic-wall case, we have instead of (2#32j) 



(2.33) 


Further, we have the condition 

v(l) = 0 . (2.34) 

If we neglect the slip velocity and temperature jump, the 
appropriate boundary conditions are, 

5(1) = 0 . (2.35) 

f 

T(1) = ^ or 5(1) = 0. 


(2.36) 
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and v( 1) = 0. (2,37) 

as before. 

There are essentially five boundary conditions in 
the free stream given by equations (2.24) to (2*28). The 
conditions (2.29) and (2.30) simply restate for r *♦ 00 the 
information contained in the algebraic equations (2.20) 
and (2. 21 ). 

On the body we have four types of boundary condi- 
tions, namely, 

(1 ) body temperature is prescribed and slip boundary 
conditions are used, 

(2) body is adiabatic and slip boundary conditions are 

us ed , 

(3) body temperature is prescribed and slip boundary 
conditions are ignored, 

and ( 4 ) body is adiabatic and slip boundary conditions are 
ignored. 

In each of these four types we have three boundary conditions, 1 
which together with the five free stream conditions define 
a consistent two-point boundary -value problem for the eighth- 
order set of differential equations we have. 

There a,re two important parameters to be prescribed, 
the -Mach number and the Reynolds number Re. Further, the 
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ratio of specific heats y , the Prandtl number Pr, and the 
viscosity-temperature law have to be given. The body tempe- 
rature it. 5 and the accommodation and reflection coefficients 
a and a should be prescribed when necessary. 


2. 6 Thin-Lay or Equations 

We reduce the basic governing equations (2.1) 
to (2.4) using the following assumptions: 

(1 ) The upstream extent of disturbance caused by the 
body is thin compared to the body radius. 


(2) Radial gradient of a given variable is large 
compared to its transverse gradient. 

However, the second assumption is used only in the viscous, 
and heat conduction terms, but not in the convection terms. 
This is done to make the equations valid even in the moderated 
high Reynolds number range. The resulting equations are, 


f \ 2 P v ^ pu ^6 pu 

(p v) + + *— — 

v, r- V / i i 


-B 


-B 


0 r- 


0 . 


B 


(2.38) 


Pr + PWr+ i|l.eZ = (i pv) 

r r r B r p •> r r 




( 






-B 
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P, 


u 


+ P T U r + P — (U 0 + V) = ( U Uj.) r , 
B a 


2 2 
and pv(h + ^) r = [ j v C^| h + "% ) r ,3 


2 r - 


( 2 . 40 ) 

(2.41) 


It may be noted that in case of Pr = 3/4? 
(2.41) has an exact integral, known as the Becker’s 


equation 
int egral 


0 1 + = constant. (2.42)) 

31 

Phis integral requires that at r = r-g, where v = 0, 
should he equal to zero. Thus, (2.42) is the solution of the 
energy equation for the case of an adiabatic wall only. In 
particular, if the wall is not adiabatic and Pr is not equal 
to 3/4, the integral given by (2.42) is not a solution* 

How , in the equations (2.38) to (2*41) we introduc 
the similarity postulates as given by the equations (2.13) 
and the nondimens ionalisat ion as given by the equations 
(2.14) to obtain the following equations: 


(p v) + 2 p (u + v) = 0, (2.43) 

P* + P v v* = (v v*) ' + if u, (2.44) 

p >2 + P -P u (u + v), (2.45) 



6 7 


(2.46) 

(2.47) 

The set of equations {2,43) to (2*47) together 
with the equation of state (2.20) and the viscosity-temperature 
relation (2.21) (to ho proscribed) , will from now onwards be 
referred to as the ’thin-layer equations 1 • The analysis of 
boundary conditions for these equations remain the same as 
for the full equations, as given in Section 2*5. These thin- 
lay or equations are first obtained by levinsky and loshihara 
(1962) who integrate them for some representative Reynolds 
numbers neglecting the slip boundary conditions. 

Our principal aim is to solve the full Kavier-Stckes 
equations for extensive ranges of the essential parameters* 

The thin-layer equations are solved only oomplementarily, 
for the purpose of comparison and finding the limit of valid- 
ity of the thin-layer assumption at low Reynolds numbers'.! 


2p 0 +p v u 


t ... (5* u ) ' 

+ p U (U + V) = — , 


and v (AT + v v ) = ^ [ y (|Zi AT + v v*) ] . 
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CHAPTER 3 

METHOD OP SOLUTION 

3.1 Introduction 

In getting extensive numerical solutions of the 
nonlinear coupled, two-point boundary -value problem formulated 
in Chapter 2, the choice of a suitable and efficient numeri- 
cal technique is of extreme importance. This senses enormous, 
computer time and additional human efforts. There are a numbe] 
of methods for the solution of two-point boundary value probl 
One of the more common of them is the so called ballistic 
or shooting integration scheme ( Lev insky and Yoshihara, 

1962, and Kao, 1964b). The shooting method involves making a 
judicious guess of the unknown initial conditions at one bouu 
dary, and performing the integration of the equations using 
these initial conditions upto the other boundary. If the pre 
scribed conditions at the other end are not satisfied, the 
initial guesses are corrected by some iterative scheme 
(for example, the method of influence coefficients, Shih anc 
Krupp, 1967 ). The task of making a suitable initial guess 
is of prohibitive difficulty in a complex system of equatio; 
as in the present case because of the sensitivity of the so 
tions to the assumed values at the first boundary. In addit 
the shooting method becomes difficult to apply when the bot 
dary conditions to be matched are specified at r 00 , as J 
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the pres cm case. 

Another method distinct from the shooting one 
is the finite-difference method. A scheme of finite differenes 
suggests itself as one means of satisfying the houndary 
conditions for a two-point bound ary- value problem. Essenti- 
ally, the application of a finite-difference scheme to the 
differential equations leads to a system of nonlinear algeb- 
raic equations, which have to be solved by some iterative 
scheme. Since the conditions at both the boundaries are 
satisfied in every iteration, the main problem in applying 
the shooting method is eliminated. However, the problem 
remains to choose an iteration procedure for the solution 
of nonlinear algebraic equations. One usual scheme is to 
apply an approximate linearisation of the nonlinear equations 
(Holt, 1964 ), and to obtain the solution of the linearised 
set by some standard method. A more direct approach is an 
iterative solution of the nonlinear algebraic equations 
themselves . Chung et al* ( 1968 ) have used the Newton's metho 
for this purpose. This method requires a very accurate initi 
guess of the variables in the entire range of integration 
in order to obtain the solution, later Dellinger ( 197 . 1 ) soli 
the same set of equations as Chung by a remarkably simple ' 
and rapid finite-difference method known as the accelerated 
successive replacement scheme. 
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The method of accelerated successive, replacement 
is originally deviced by Lieberstein (1968) for the iterative 
solution of 'mildly nonlinear* elliptic partial differential 
equations. This is a generalisation of Young 1 s (1954) method 
of over-relaxation for the iterative solution of the linear 
algebraic systems derived from linear elliptic partial diff- 
erential equations. This method is first applied to the 
boundary-layer equations by Lew (1968) who solves the Blasius 
equation and Strom (1968) who studies three-dimensional boun- 
dary layers, later Dellinger (1971) applies this method in 
a modified form to nonequilibrium stagnation flows under 
hypersonic conditions. There are two salient features of 
this method. One is that it involves the application of one 
particular equation only for the correction of one particula: 
unknown variable in each iteration. The second is that the 
corrections applied to the values of the variables at each 
of the mesh points are controlled by 'acceleration factors', 
which prevent the iteration scheme from diverging. Thus, a 
successful application of this method does not critically 
depend on how well the initial guesses approximate the con- 
verged solution. Because of these advantages of the acclera,- 
ted successive replacement scheme, we find it the most suit 
able for the present investigation. J 
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3.2 Accelere.t ed Successive Replacement Scheme 

The equations to be solved, either the full Navier- 
Stokcs equations (2.15) to (2.21) or the thin-layer equations 
(2.43) to (2.47), are first expressed in terns of a trans- 
formed normal coordinate n which is defined by 


n = 



(3.1) 


where n = r - 1 and n Q = r e - 1. r Q represents the effective 
nondimens ional free stream radius, that is, the radius at 
which the flow variables have virtually -attained their 
proper asymptotic free stream values, r and hence n are 
apriori unknown, but are determined as an integral part of 
the solution. The transformation (3.1) keeps the body at 
n = 0 and the free stream at n = 1 . 

Having made the above transformation, one observes 
that the unknown constant n Q occurs explicitly in the trans- 
formed differential equations. The technique of accelerated 
successive replacement is applied only for the second-order 
equations, while the first-order equations can be solved by 


direct numerical quadrature. There are three different ial } 
equations which are of the second order, which may be writte 


as 


U (u ', u\ y , u, v , p , p, u, p 9 , n ) = 0, 


f 


(3.2) 
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-n _ r _ _ t _ _ __ 

^ (v , u , y’, p , V, p , y , nj = 0 , (3.3) 

T (T » T , y» , p', u, v, p* , v , n e ) = 0, (3.4) 


prime denoting differentiation with respect to n . Solutions 
are desired in the range 0 < n <_ 1, which is divided into 

-t 

iTDIY number of equal intervals, each of length A = 

Equations (3.2) to (3.4) are replaced by their f inite-diff er 
once approximations . Eor the first and second derivatives 
we use, 


% 


% f 1 ~ %-1 

2A 


(3.5) 


tf 

0vXl cl ~ 


X N +1 ” 2 % + %-1 

~ — 2 

A 


(3.6) 


respectively, where stands for the value of any one of 
the variables u, v, T, v , p at the Nth mesh point. N = 1 
corresponds to the body, and N = NDIY +1 to the free stream 
edge. To start with, initial guesses are made for all the 
flow variables at each of the mesh points, consistent ■with, 
the boundary conditions. To correct the variables u, v and 
T, each is associated with a particular equation from the 


set of equations (3.2) to (3.4). To generate corrections .to 
a given variable only the corresponding equation is used. 

It is natural to associate u with (3.2), v with (3*3) and 
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T with (3.4), where these variables respectively occur in 
their second derivatives. 

We illustrate here how the corrections to the 
initial guesses are uie.de for the variable v. The other vari- 
ables u and T are similarly corrected. The (k+1) approxi- 

— ijii 

nation to v at the N mesh point is calculated from 






(p|L) 

V dv 1 


N 


(3.7) 


The number o) v , called an acceleration factor, is calculated 
for each correction in the following way. First, it is 
required tha.t 


v- 


;k+1 

N 


- v 


-k 

N 


^N 


v 


Y 


N 


^ av. 


I 


dv w 


< e 


(3.8) 


where 6 is a prescribed, snail, positive number. Now, we 
define 




e 

Y 


N 


-k 

t n 


dY 


N 


dv 


N 


(3.9) 


"Whenever a < 1, we put and a v > 1 we put to y = 1 • 

Thus, the acceleration factor ensures that the correction 
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applied to a variable is never greater in magnitude than 
6-times the previous value of the variable, The correction 
starts from the second mesh point, that is, N = 2, to the 
penultimate point, that is, N = NDIV. Once a variable is 
corrected its previous value is erased, and the new value 
is used in any subsequent calculation. 

Having corrected u, v and T once, at all the mesh 
points between the body and the free stream, the remaining, 
variables are calculated as follows. The viscosity yT is cal- 
culated from the new values T. p and p£ are calculated by 
numerical quadrature from the free stream side to the body, 
using the equations in which p" and pg occur in first deriva- 
tives. Pressure is usually determined from the equation of 
state. This completes one cycle. The next iteration is taken 
up by once again correcting u, v and T by accelerated succ- 
essive replacement and then p" and pg by quadrature and so on, 

The general effect of this procedure is to decouple 
the system of equations in each single correction. The effect 
of coupling of the equations comes through the successive 
iterations that are carried out* 

To perform even the first iteration, the effective 
free stream location, n Q , should be known. To start with, 
some guess for the value of n 0 is made. Provided a proper 
value of n is used, the procedure given above can be iterated 
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to obtain the desired solution. However, if the value specified 
for the free stream location is too small, the solution does 
not show proper asymptotic behaviour at the free stream, the 
profiles there having appreciable first derivatives. On the 
other hand, if the value assumed for n 0 is larger than actually 
needed, difficulties can arise in maintaining a stable calcula- 
tion in the uniform free-stream region, at least in the high 
Reynolds number case. Moreover, if a large portion of the inter- 
val (0,1) of n is to be occupied by a uniform free stream, 
then the effective calculations are done only in a lesser number 
of mesh points, thus lessening the accuracy of the calculations. 
These situations are avoided by continuously correcting n in 
each of the iterations, or once after a certain prescribed 
number of iterations. The criterion to correct n g is a proper 
asymptotic behaviour of the variables. Since temperature is 
found to be the most receptive of the impending disturbance, 
as one comes from the free stream towards the body, the value 
of n 0 can be adjusted to a requirement that the first derivative 
of temperature at the free stream should be small. 

3 • 3 Full lavier-Stokes Equations in Finite-Difference Form 

After making the transformation (5.1) on the full' 
i'Tavier-Stokes equations (2.15) to (2.19) we obtain the following 
set of equations i 
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(1 + n e n )v 


(1 + n n) 

{~ 2 u - 2 v =— 2 . y } , ( 5 . 10 ) 


8 E + v £ + ^ / u + 1 ) 

^ (1 + n e n ) n e ^ 1 + n e ri 


Re p u 


(■AVt ) + 


+ ; 1 + » e ” 


Re p u y 
n e * 


- 1 

v _ 

3n e (T”+ n e n. 


il ( 2. + _1R) + i22E!_ 

n e H + n e n n e (l + n) - 


(3.11) 


JL .. _ 3 Re (E 

nf“ “ 4 i 


! < - t -» 

+ P v v\ / 2 , y \ jv 

^ y — \,T' 'll w ’ r n — ' 


n y v 1+nn n 7T y n 

e e 1 e ^ e 


u , r 7 , ^ K } ( . u . * v 

2 n e (1 + n e V) 2(1 + n n ) . ■ n y 1 + n e n 

0 © 


(3.12) 


EL— + 


TT+l i e n 


(u + v). 


(3.13) 


and 
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— — 2 — t — 2 — _ ? 2 2 

A p v (1 + n p n) T v P (l + n p n) 2 v v (l + n_ n ) 

^ 


a u - _ 2 p iT (l + n n )v * 2 

T5e“ (« + ^) - Ji£ { 3T + 2 (u + v) > 

© 


A(1 + n p n) _ — 1 — 1 — j — — « 

+ — — — ~p— — - {2n p I + (1 + n n)( v I + VT)} 

Re Pr n c . e 


(3.14) 


It may be noted that a prime denotes differentiation with, 
respect to n in the equations (3.10) to (3.14), while it 
denotes differentiation with respect to r in Chapter 2. Consis- 
tnetly, throughout Chapter 3, a prime denotes differentiation 
with respect to n . 

The eolations (3*11), (3.12) and (3.14) which contain 

second derivatives of the variables are now put in a finite- 

"fh 

difference form at the N mesh point using the difference 
formulae (3.5) and (3.6). We obtain, 

TT %+1 “ 2u h + U M-1 8 U N + 

j-r + y FT^FTTF 2 

0;: + i " wh— i ) h, + F 

+ 24 n„^ 1 + n e (N-1)A + 
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Re p N U I'T r U R + V N 


u t 


f u 1 

1 1 + 3“ (1-1 ) A J + 


Re P E ~ 

n e 2 1 " 
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( V I + 1 " v H-i) 


6 A n e {1 + n e (N-1)A> 


(% + 1 u N-i) 
2 A n “ 


1 +“n e (N-l)A 


( y R+1 ~ 

2A n. % 


+ 


2 p 2K Re 


\1 + n p (17-1 ) A\y 


= 0, 


(3.15) 
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t H+ 1 - 2v N + V N-1 
~ 2 2 
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5 Re { 1 “ PN ” 1 1 " V N- 1 ^ } 
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2 A n 


1 + n (11-1 )T + 2 A n u.r 2A n 


(Pl\T , -i - Put ) , ( v 7 
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e y N Ai e 
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4 A i n Q (R-1 ) A+l > t n e 
6 
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= 0, 


(3.16) 
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2 A n 


A P N V N {1+n e ^ N ~~' ! ^ A} ( T N+1 ~ T N-1 ^ Q (3.17) 

~ ^ ~~ Hfi ‘5 A 


Prom the above equations one obtains the derivatives 


dU 2 8 U N+1 “ y N-1 

d % A 2 n 2 ” 3{1 + n e (S-T ) a> 2a {1 + n e (N-l)A} 

Re (2 u F + v N ) 

wit {1+n e (N-l)A> 


f 


(3.18) 



80 


av = 2 

dTj, i 2 n 2 


3Re p h - v^ ) 

8 4 n e ^ 


+ { 



1+n (IT— 1 ) a 


%+1 ~ U N -1 } 
2L n e 


■ 1 

{1+n (R-l) a } , 


and 



2 

2A u N {l+n e (N-l) A } 

Re Pr n 2"*^ ” 

e 


(3.19) 

( 3 . 20 ) 


These derivatives are used in determining the acceleration 
factors ^ U 9 w o)^ at each step, 

3 • 4 Thin-Layer Equa t ions in Finite Difference Porm 

The thin-layer equations (2.43) to (2.47) give* after 
the transformation (3.1 ) 


(p v) + 2 n P (u+v) = 0, 


-< 2 -- 


2 ng p 2 + n e P v u + n g pu (u+v) = 
4 

3Pe 


(u u) , 


(3.21) 

( 3 . 22 ) 


^ — (v v*) = n 0 (p T +p v v* 


2 - - f 

3 Re 


v U + ^Ra u U )»(3.23) 


- r _ t 

*2 


p 9 + p - n e 


p U (u + v) = 0; 


(3.24) 
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and 

pv n e (A T 


_ _ t 

+ V V 


+ ""’ } ’ 


(3.25) 


form, 


Keeping tie second-order equations m 
we obtain, for the H th mesh point, 


f init e-dif f erenc e 


_ y, - 2 

% = Ke~ V H ~ k 2 

Js±iJJd=i_ 

+ ^5 2Z 2 ‘ 

- - Vi ~ Ult ~ 1 

- 2 n e P 2rT - n e P H 2 t 

fc. t-Jf 

_ (3.26) 
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N 3 Re 
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3 Re 


2A 


2A 


- - n„ . - - T B+1 ' ‘ 

^jti - 3 =n. p b v H 24 

e ' ^ 


- (—”9 A 


O — 11+ 1 

rne ‘ 


2A 


V,- • “ U ’ 

+ 3~Re 2A 


= 0, 


(3.27) 
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3.5 Boundary Conditions in Finite Difference Form 

a) In all the cases the free stream boundary conditions 
are as follows: 


U NDIY+1 ~ 1 

V HDIY+1 = _1 

- 

%DIV+1 = 1 ~ 2A 

P p = 0 

rlDIV+l 


(3.32) 


P NDIV+1 ” 1 

b) On the surface of the body, that is, at the first 
mesh point (denoted by the suffix 1), the slip boundary condi- 
tions are as follows : 



(3.33) 


and 


T = T + 
1 x w + 


2y 2-a 
T +1 a 



^ dl \ 

Pi/ T 1 Re ld?7 r=1 5 


(3.34) 


where the derivatives on the rigbt-hand side of these two 
equations can be expressed as 
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r=1 


4 u 2 - u 5 - 3 Ul 
2A n e 


and 


I dT 1 


\dr / 


r=1 


4 T 2 - - 3 T 1 

2A n e "" 


(3.35) 


If the slip velocity and temperature jump are 
neglected, we have simply 


= 0 


T„ = T 
1 w 


(3.36) 

(3.37) 


In the adiabatic-wall case, instead of the condition 
(3.34) or (3.37), we have 


T 


4 ^ 2 ~ 


(3.58) 


1 3, 

In all the cases, the normal velocity vanishes at 
the body, which means 


v 1 = 0 _ (3.39) 

While the conditions at the (HD 17+1) mesh are 
fixed, the conditions at the first mesh point are in general 
not fixed, and the values of u-j and T 1 are to be changed 
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after each, iteration. 


3 . 6 Quadrature for p and P o 

After each round of applying accelerated successive 
replacement to the values of u, v and T, new values of p and P2 
are obtained by numerical integration of the first-order 
equations (3.10) and (3.13) or (3.21) and (3.24) respectively 
depending on whether we are solving the full Mavier- Stokes 
equations or the thin-layer equations. Since the values of P 


and Pg are known at the free stream* but unknown at the body, 
the integrations are carried out from the free stream side to 
the body surface. The values of p* and P£ at the first point 
before the free stream end are calculated using the trapezoidal 
rule. The values at the next four mesh points are calculated 
using the Simpson's rule. At all the remaining points upto the 
body surface the Weddle’s rule is used. 

In case of the full Navier-Stokes equations, the 
integrals to be evaluated are as follows: 


r K r 

/ -rr dn = f 


n 


(1+n 0 Ti )v 


0 _ (1+iyO _ r 

2u - 2v - e 


n„ 


v 


} dn 


(3.40) 
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and 

/ p~2 an-n e / t eLf an (3.41) 

6 G 

Equation (3.40) presents a difficulty at the wall 
where v = 0. The integral on the right-hand side of the 
equation (3.40) cannot he directly calculated at the body^ 

which is the last mesh point for the density integration. 

p 1 

However, the value of - = — ■ at the body surface may be found 

p 

as follows. Assuming that there is some, non-infinite normal 
gradient of density at the body, which appears to be true on 
physical grounds, it follows that the numerator of the integ- 
rand on the right hand side of equation (3.40) should vanish 
at n = 0. Thus, one gets at n = 0 

v 1 = - 2 n e u 1 (3.42) 

* 

Hence v^ is non-zero in the case of slip boundary conditions 

and the value of Ll at the body may be evaluated after applying 

p 

the L 'Hospital rule to the integrand. In case of no-slip 
boundary conditions, an application of the I’Hospital rule twice 
is necessary. However, it is not essential to use the L'~ 
Hospital rule at all, because, using the equation (3.12) and 
the differentiated form of (3.10) one can easily obtain an 
explicit expression for p at n = 0, which can be used to 
integrate p upto the body surface. 
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In the case of a very cold wall, where the density 
is expected to have a very large gradient at the body surface* 
integrating the equation (3.40) upto the body Is not very 
convenient. This is because* the value obtained for p" at the 
body* is used to calculate the pressure on the body, in in- 
accurate pressure* calculated from an inaccurate density, 
affects the subsequent calculation of the corrections to the 
normal velocity v* which in turn greatly affects the density 
near the wall. A more convenient method to predict the wall 
pressure, in the case of a very cold wall, is to integrate 
(3.40) only upto one mesh point away from the wall. The pres- 
sures upto this part are calculated from the equation of state. 
The final pressure on the wall is obtained from the linear 
expression 

- _ t 

P 1 = P 2 “ & P 1 (3.43) 

_ t 

where the derivative p^ can be obtained from writing the 
equation (3.12) at the first mesh point, using appropriate 
finite-difference approximations to the derivatives appearing 
in that equation. This approach avoids difficulties associated 
Yv'ith integrating (3.40) in a region where the density is varying 
rapidly. Also, the use of an 6 - criterion, such as in the 
equation (3.8), for the maximum relative change in the value 
of pat any particular station, Is found to be of great use 
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in controlling the instability and quickening the convergence. 

3 • 7 Convergence of the Num e rical Procedure 

The numerical procedure described in the previous 
sections is programmed in the Fortran IV language and the 
computations are carried out on an IBM 7044 computer at the 
Indian Institute of Technology, Kanpur. Table 3.1 indicates 
the approximate computer time required for various operations. 

Table 3.1: Approximate Computer-Time Required for 
Various Operatio ns 


8 .No. Operation Time in Seconds 


1 . Compilation of the programme to 140 

solve the full Navier-Stokes equations 

2. Compilation of the programme to solve 120 
the thin- layer equations 

3. Execution of 500 iterations to solve 340 

the full Navier-Stokes equations 

4. Execution of 500 iterations to solve 250 

the thin-layer equations 

5 . Execution of 500 iteration to solve 200 

the thin- layer equations in the 
insulated-wall case using Becker's 
integral (Eqn. (2.42) ) . 


In general, the interval 0 <n < 1 is divided into 
60 equal sub-intervals. The maximum relative correction 6 in 
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equation (3.8) is taken to "be equal to 0.01. 

Figures 3.1 and 3.2 show the convergence of the 
pressure and normal velocity profiles respectively, in a typi- 
cal low Reynolds number case. The initial profiles are arbi- 
tarily prescribed, but they are consistent with the boundary 
conditions. In the first fifty iterations, the acceleration 
factors a> u j as v , airj, (equation (3.8)), are less than unity at 
most of the stations. This means that the corrections are still 
being controlled. By about 100 iterations, all the acceleration 
factors are equal to unity, showing that the relative corrections 
are by themselves less than £ = 0.01. By about 500 iterations, 
the profiles converge to a, limit. At this stage the maximum 
correction applied to any variable has an effect only in the 
fifth significant digit. Calculations carried upto 1000 itera- 
tions do not show any graphical difference from the profiles 
obtained at 500 iterations. 

Regarding the initial profiles to be prescribed, it 
is enough if they are consistent with some of the boundary 
conditions and do not have abrupt jumps at least near the body 
surface. The numerical procedure has the capacity to assimilate 
even large initial discontinuities on the free stream side* 

This is shown in Figure 3.3, where a moderately high Reynolds 
number case is studied. The initially prescribed temperature 
profile is not consistent with the free stream boundary 
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condition for the temperature. This considerable initial 
jump is smoothened out in course of iterations, and within 
300 iterations we get a converged profile. It may be noted 
that, in the case presented in Figure 3.3? the effective 
radius r e (equation (3.4)) is kept fixed at 1,4. In the final 
converged solution, this results in a uniform free stream 
zone of about 40 percent of the total range of calculations. 
Thus, of the total number of 60 points of division, only 
about 35 correspond to the effective profile. It is verified 
that a decrease of the effective radius from 1,4 to 1.25 does 
not alter the converged solution. 

We have also verified, by performing the integrations 
for some individual cases, that 

(i) increasing the number of divisions of the interval 
0 < n< 1 from 60 to 120, and 

(ii) increasing the maximum allowable relative correc- 
tion 6 from .01 to .05 do not have any appreciable effect on 

.0 

the converged profile. - 

Once a converged solution is obtained, for a parti- 
cular set of the parameters like and Re, this solution can 
be prescribed as the initial guess for adjacent values of the 
parameters. Usually, with a 20-50 percent change in the 
Reynolds number, one can obtain a new converged solution in 
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about 400 iterations. Profiles are supposed to have converged 
if for each further iteration variables change almost in a 
fifth significant digit. 

As discussed so far, the coefficients to , and 
which are calculated at each mesh point in each iteration, 
are used only to control the corrections applied (according 
to the condition (3.8)). Phis control of corrections is very 
important in the initial iterations. If, instead of using the 
condition (3.8), we always use unit values of w , and uj 

the iteration procedure may diverge even in the initial . stages. 
However, it is not necessary to impose the condition (3.8) in 
all the iterations. This is because of the fact that by about 
100 iterations the profiles would have adjusted themselves 
well enough to the given parameters and the condition (3.8) 
is by itself satisfied even without resorting to the accelera- 
tion coefficients , oj . In such a case, we can further accele- 
rate the convergence by taking values of greater than unity. 
Por example, we can impose the condition that u = 1.5 if 
'L 1 *5 (cf. equation (3.9)). This method is tried in some 
particular cases, and it is found that there is a considerable 
gain in the number of iterations required to achieve a given 
measure of convergence. 
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FIG. 3.2 -CONVERGENCE OF NORMAL 
VELOCITY PROFILE 


FIG. 3.3 -CONVERGENCE OF THE TEMPERATURE 
PROFILE, SMOOTHENING OF AN INITIAL 
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CHAPTER 4 

RBSUITS FOR THE COLD-BODY CASE 


4 * 1 Introduction 

In this chapter we discuss a wide range of the cold- 
body solutions, spreading from the boundary- layer limit to 
the transitional -flow regime. In particular; 

(1) The effect of slip velocity and temperature jump 
on the detailed flow field and the overall characteristics 

of the flow is investigated, 

(2) The validity of the thin-layer assumption is 
discussed, 

(3) Comparisons are made with the results based on 
the theoretical work of Cheng (1961, 1963a, 1966), levinsky 
and Yoshihara (1962) and Liu (1967), 

( 4 ) The solutions are compared with the Monte Carlo 
calculations of Vogenitz and Takata (1970), 

(5) Detailed low- Reynold s-number flow studies are 
made in order to know the effect of wall temperature and Mach 
number variation, 

and, (6) Extensive comparisons are made with the experi- 
mental results for the stagnation-point heat-transf er rates, 
and the detailed density and temperature profiles* 
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Humerical solutions are obtained for the following 
sets of parameters? 

I * (a) y = 5/3, *L= 10, Pr = 3/4 

y = °- 029 V 

( c ) 1 < Re < 300, in close steps of Re 

(d) va fl 

(e) a = a = 1 , 

and 

II. 

(a) Y = 1.33, = 10, Pr = 0.71 

0>) T w - 0.1 Iqoq 

(e) 1 < Re < 300, in close steps of Re 

(d) p a fl 

( e ) a = a=1. 

Further, 

III. Retailed case-studies of the effect of Mach number 
and the wall temperature variation are made at a typically 
low Reynolds number. These solutions are discussed in 
Sections 4.5 and 4.6, 
and 

XV. A number of individual cases are solved for the 
purpose of comparing the detailed flow profiles with experl 
ments. These comparisons are made in Section 4.7. 



97 


4.2 General Plow-Profiles 

All the solutions presented in this section are 
obtained using the full Navier-Stokes equations (2*155 to 
(2*21) with the slip boundary conditions, Figures 4*1 to 4*5 
are from the set X. Figures 4*6 to 4.10 correspond to the 
set II* 

In Figure 4.1 , we draw the profiles for the normal 
and tangential -velocities, the temperature T and the pressure 
p, at Re = 296, The extent of disturbance caused by the body 
is about one fifth of the body radius. The normal velocity, 
temperature and pressure have steep variations in the region 
1,17 < r < 1*19* This corresponds to the shock wave. Hear 
the body surface, r = 1 , the temperature and the tangential 
velocity u have steep gradients and show a boundary-layer 
nature. From r = 1.0 to r = 1.04, which corresponds to only 
one fifth of the total extent of disturbance, the tangential 
velocity attains about 80 percent of its free-stream value 
and the temperature attains about 95 percent of its maximum 
value. Further, the variation of pressure is negligible in 
the region 1.0 < r < 1.1. The flow In the region 1.06 < r 1.16 
can be considered as essentially inviscid. The slip velocity 
and the temperature jump are very small. 

The profiles at Re = 34 are drawn in Figure 4*2. At 
this Reynolds number, the extent of disturbance considerably 
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increases compared to the previous case, and is about one 
third of the body radius* We do not find abrupt variation 
in the flow quantities in the shock-^wave-like region. The 
viscous layer near the body and the shock-wave-like region 
merge into each other and the inviscid region separating them 
no longer exists* The pressure has a considerable gradient 
near the body surface, and therefore, the usual boundary-layer 
approximations fail in this case. The maximum $ in the flow 
is about 0.93 (T/T^ =32) against 0.97 in the previous high 
Reynolds number case. This decrease in temperature is due to 
the weakening of the shock wave. Consequent to merging the 
effect of cooling the wall is propagated to the shock-wave 
like region as well. The inflexion in the density, p , near 
r=1.15 is due to the shock compression, while the steep incr- 
ease in near the body surface is due to the strong cooling 
of the body. 

In Figure 4.3, the profiles at Re=3*4 are presented. 

At this highly rarefied flow condition, the extent of the 
disturbance caused by the body is even greater than the body 
radius. One cannot expect any type of a thin-shock-layer 
approximation to be valid to describe such a flow. Here, the 
pressure profile is much different from the high Reynolds 
number profiles in that the gradient is large near the body. 
The shock-wave-like region thickens appreciably and is indis- 
tinguishable from the viscous flow near the body surface. The 
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slip velocity is almost 60 percent of the free— stream tangen- 
tial velocity and the temperature of the gas at the wall is 
about 15 times the free-stream temperature. The overall com- 
pression of the gas is very much reduced, the density of the 
gas at the body being only 12 times the free-stream density, . 
against a multiple of 42 at Re = 34. The maximum temperature, 
T maZj attained in the flow is 0.78 {T/T to - 27.0) in the present 
case, as compared to 0.93 (T/T^ =32.0) at Re = 34 and 0.97 
(T/£ eo *33.0) at Re = 296. 

In Figures 4.4 and 4.5, we present the composites of 
the normal and tangential-velocity profiles respectively, 
at various Reynolds numbers. From the normal-velocity profiles, 
the rapid thickening of the shock-wave- like region, as the 
Reynolds number decreases, can he clearly visualised. Very 

near the body the normal-velocity profiles do not change 
appreciably with Re, compared to the change in the shock- 
wave-like region. However, the tangential-velocity profiles 
show a considerable difference, indicating the thickening 
of the viscous layer near the body. F 0 r a given Re, the tan- 
gential velocity attains its free-stream value at a smaller 
value of r than the normal velocity. This is so because, in 
the shock-transition zone the tangential velocity does not 
suffer as much change as the normal velocity. 
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Figures 4.6 to 4.10 are drawn from the solutions belo- 
nging to the set II. The basic features of these profiles are 
same as those of set I. Ik. 

i, ' *. v • t 

In Figures 4.6 and 4.7 we draw the flow profiles for 
Re = 191 and 3.6 respectively* At Re = 191, we observe a 
moderately thick shock, but still separated from the viscous 
layer near the body by an essentially inviscid region* The 
inviscid region corresponds to the plateau in the density and 
temperature profiles. Coming from the free-stream side, the 
initial rise in density is due to the shock compression. After 
having almost a constant value in the inviscid region, the 
density again increases rapidly due to body cooling. The 
variation of pressure near the body surface is snail, showing 
that a classical boundary-layer type of flow exists near the 
body. 

In Figure 4.7, the fully-merged profiles at Re =3*6 
are drawn. Pressure has a steep variation near the body* The 
density profile (not shown in the figure) behaves In a simi- 
lar way. The slip velocity and the temperature jump are con- 
siderably large compared to the previous high Reynolds number 
case. The maximum temperature in the flow is reduced. At 
Re = 191, T ma;£ = 0,98; and at Re = 3.6,T maz = 0.63. As men- 
tioned earlier, this is due to the interaction and merging 
of the shock-wave-like region with the viscous layer near 
the body. 
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In Figure 4.8, temperature profiles are drawn at four 
different Reynolds numbers. The decrease in the maximum tem- 
perature due to merging, the thickening of the viscous layer 
near the body and the more rapid thickening of the shock-wave- 
like region are clearly visualised in this figure. Figures 
4.9 and 4.10 give the pressure and the tangential velocity 
profiles respectively, at various Reynolds numbers. As already 
mentioned, the pressure gradient in a region near the body 
increases appreciably as the Reynolds number decreases. 

Figures 4.8 and 4.10 show that the temperature jump and. the 
slip velocity increase as the Reynolds number decreases. The 
effect of slip boundary conditions on the flow is discussed 
in detail in the following section, 

4 . 3 Effect of Slip Boundary Conditions 

In this section, we investigate the effect of the 
slip velocity and temperature jump, by comparing. the solutions 
obtained using the slip conditions (2,31) and (2.32) with 
those obtained using the no-slip conditions (2,35) and (2,36). 
Throughout this section we present the solutions of only the 
fully Navier-Stokes equations. 

In Figure 4.11, we draw the velocity and temperature 
profiles at Re = 170, for the case of Y = 1.667, = 10 and 

T w = T to (T w = .029). At this Reynolds number, except a slight 
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difference near the body, the slip and no-slip solutions 
are almost identical in the entire layer of disturbance. 

In the slip case, the predicted tangential velocity at the 
wall, u(1), is 0,0118. The temperature of the gas near the 
body is 1,67 times the temperature of the body. This rise 
in temperature is small compared to the high level of tempe- 
rature that the gas attains even at a small distance from the 

body. 

Figure 4.12 presents an entirely different situation. 
Here, the velocity, pressure and temperature profiles are 
drawn at a typically low Reynolds number, Re = 3.4, while 
the other parameters are same as for Figure 4.11. It is. obser- 
ved that the slip boundary conditions change the entire 
structure of the flow in the shock layer upto the free stream. 
The temperature predicted by using the slip conditions is 
always greater than the no-slip value. The pressure level near 
the body is decreased with slip, while at some distance from 
the body the pressure increases over the no-slip solution. 
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Incorporating the slip conditions slightly increases the 
extent of disturbance caused by the body. F or example, the 
temperature almost attains its free stream value at r = 2*15 
in the case of no-slip solution, while with slip boundary 
conditions, this is not so until r = 2.4. The slip velocity 
and the temperature jump at the body are large. Slip velocity 
is about 28 percent of the corresponding free stream value. 

The temperature of the gas at the wall, 5(1 ), is equal to 
0.438, compared to T w of 0.029, 

Figures 4.13 and 4.14 are from the set II of computa- 
tions. The observations made from Figures 4.11 and 4.12 are 
valid for this case as well. At Re = 71*7, shown in Figure 4.13, 
the effect of slip boundary conditions is confined to a small 
region near the body. But, at Re = 7*17, as shown in Figure 
4.14, the slip conditions affect the profiles even upto the 
free stream. At Re = 71.7, the slip velocity is 0.035, while 
at Re = 7.17, the slip velocity is 0,152. The temperature, 

5(1) of the gas near the body is 0.151 at Re = 71.7 and 
0.297 at Re = 7.17. Thus, with a tenfold decrease in the 
Reynolds number, the slip velocity and the ratio of tempera- 
ture jump to the wall temperature are increased only by mul- 
tiple of four. However, it is important to note that at low 
Reynolds numbers, the entire shock layer is affected appreci- 
ably, while at high Reynolds numbers this effect is small 
and is confined to a small region near the body. 



104 


In Figure 4.15, we draw the heat-transfer rate to 
the hod 3 ?' versus Reynolds number. The heat-transfer rate 

0 rjnt i 3 

(k 9 ^-) w is nondimensionalised with p °o V, * Slip and no- 
slip curves are drawn for the sets I and II. In both the cases 
it is observed that the effect of slip is to increase the 
heat-transfer rate to the body. When the slip boundary condi- 
tions are used, there are two. opposing effects. The tempera- 
ture gradient at the wall is decreased compared to the no- 
slip case. But, the temperature itself is increased, thus 
increasing the coefficient of thermal conductivity, k. The to- 
tal effect happens to increase the heat-transfer rate at a 
given Reynolds number Re. In the case of y« 1.33, this incre- 
ase is small even in the very low Reynolds number regime. 

But, in the case of y = 1.667 the increase in heat-transfer 
is appreciable at low Reynolds numbers. As the Reynolds number 
is decreased the difference between the slip and the no-slip 
solutions tends to increase. 

Figure 4.16 shows the variation of the skin-friction 
coefficient (^f|) w / p oo • For Y = 1.667, the effect of slip 
conditions is considerable. For y = 1.33 this effect is compa- 
ratively very small. In both the cases slip increases the 
skin friction. As the Reynolds number decreases the difference 
between the slip and the no-slip cases tends to increase. 

In Figure 4.17, we present the variation of the slip 
velocity and temperature jump with Reynolds number. Curves are 
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drawn for the sets I and II. In both the cases, the slip 
velocity and temperature jump increase rapidly as the Reynolds 
number decreases. In the higJi-Reynolds-number limit they tend 
to the corresponding no-slip values asymptotically* Ror Re 
greater than about 25, the slip velocity for r= 1*667 is 
only slightly less than that fory = 1.33. But in the low- 
Reynolds-number regime, the monatomic gas shows an apprecia- 
bly larger slip velocity than for the case of y = 1.33. Ror 
example, at Re = 5, with y = 1.667 the slip velocity is 0.240, 

while with y= 1.33 the slip velocity is only 0.174. This 

change in the slip velocity is due to the change in the value 
of y, rather than to the difference in the wall temperature. 
We also remark that 3(1)/^ is. larger for the case of 
1 = 1 .667 than forY= 1.33 for all Reynolds numbers. 

In epitome, we make the following comments regarding 
the effect of slip boundary conditions on the flow past a 
cold body: 

(1) At low Reynolds numbers, the slip effects are 
large and the slip boundary conditions affect the entire 
flow field including the shock-wave-like region. 

(2) The extent of disturbance is slightly increased 




due to slip 
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(3) Slip boundary conditions increase the skin 
friction and the heat-transfer rate to the body, 

(4) The effect of slip conditions depend intrinsi- 
cally on the ratio of specific heats Y , lor a larger y the 
slip effect is larger. 

and, 

(5) The effect of slip conditions decrease as Re 
increases. It is negligibly small at Re = 300. 

4.4 Limit of Validity of the Thin-layer Assumption 

At low Reynolds numbers the extent of disturbance 

caused by the presence of the body is of the order of the 

body radius, and the assumption that r* r^ in the layer of distur- 
■valid 

bonoe-isnot . .The solutions obtained from the thin-layer 

A 

equations (2.43) to (2.47) are compared with the solutions 

of the full Navier-Stokes equations (2.15) to (2.21). Compa- 
risons include slip as well as no-slip cases. All the results 
presented in this section are for the parameters Y = 5/3, 

Moo = 10, Pr = .75 and T w = T^. 

In Figure 4.18, we compare the temperature profiles 
calculated from the full Navier-Stokes equations and the thin- 
layer equations, at five different Reynolds numbers. Re = 6.8, 
10.2, 17.0, 34.0 and 85.0. All these profiles are for the no- 
slip case. At Re * 85.0, the full Ravier-Stokes and. thin-layer 
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solutions give exactly the same results except for a small 
difference near the peak temperature. The maximum temperature 
predicted by the thin-layer equations is about one percent 
less than that presented by the full Ravier-Stokes equations. 
At Re = 34.0 this difference increases to about 6 percent. 

But still the flow near the body and in the shock-wave-like 
region are almost the same in both the cases. At Re = 17.0, 
we observe that the profiles are different in the entire 
range except very near the wall. The thin-layer equations 
predict a larger extent of disturbance than the full Ravier- 
Stokes equations. At Re = 6.8, there is absolutely no simi- 
larity in the profiles of the thin-layer and the full Ravier- 
Stokes equations. The thickness of the thin-layer profile 
is about three times larger than the thickness of the full 
Ravier-Stokes profile. 

Figures 4*19 and 4.20 are complementary to Figure 4*18. 
In Figure 4.19, we draw the variation of the maximum temp- 
erature in the flow field with Reynolds numbers. The differ- 
ence between the thin-layer and the full Ravier-Stokes cases 
increased appreciably as the Reynolds number is decreased. 
Compared to the predictions of the full Ravier-Stokes 
equations, the thin-layer equations differ by about 2 percent 
at Re = 100, and about 4 percent at Re = 50. At Re = 20 this 
difference increases upto 10 percent. In Figure 4.20, where 
the extent of disturbance is drawn against Re, we observe 
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even more remarkable difference between the thin-layer and. 
the full Navier-Stokes cases. This extent of disturbance 
is calculated from the temperature profiles and is arbitrarily 
defined as the nondimensional distance from the body at 
which T ;= 0.1, A distance calculated from the density profiles, 
would understimate the extent of disturbance. The present 
definition of the extent of disturbance should give a reason- 
able measure of the total shock-layer thickness. Figure, 4*20 
shows that the thin- layer equations give a larger thickness 
than, the full' Navier-Stokes equations. Compared to the pre- 
diction of the full Navier-Stokes equations, the thin-layer 
equations differ by about 2 percent at Re = 50, 12 percent 
at Re = 25, 20 percent at Re = 17 and as much as about 
100 percent at Re = 8, 

In Figure 4*21 we present the pressure of the gas at 
the wall versus Reynolds number. The slip as well as the no- 
slip curves are drawn for the full Navier-Stokes equations 
and the thin-layer equations. T his figure shows. the intrin- 
sic failure of the thin-layer equations in the low-Reynolds- 
number regime. At very high Reynolds numbers all the four 
curves tend to the same value, which is the value predicted 
by the inviscid theoiy. (In the boundary-layer limit, the 
pressure variation across the boundary layer is negligible) . 

As the Reynolds number decreases, the full Navier-Stokes 
equations predict an initial decrease of the pressure at the 
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wall. As the Reynolds number is decreased further, the 
pressure increases steeply. This steep rise in the pressure 
after an initial dip is observed by Sherman (1953) and Potter 
and Bailey (1963) in their extensive impact-pressure measu- 
rements. The general nature of these experimental results 
is already shown in Figure 1,5-. In a subsequent section our 
pressure predictions are compared with the experimental data 
-and it is found that they show a good agreement both quali- 
tatively and quantitatively (see Figure 5.32), It is important 
to note that the thin-layer equations fail to show the rise 
of pressure at low Reynolds numbers, is a matter of fact. 
Figure 4.21 shews that the pressure given by the thin- layer 
equations continuously decreases, and at low Reynolds numbers 
decreases more rapidly, Levinsky and Yoshihara (1962), who 
integrate the thin-layer equations neglecting slip, conjecture 
that this continuous drop of wall pressure is due to the 
neglect of slip. Our present results show that introducing 
slip boundary conditions further decreases the wall pressure, 
be it the thin-layer equations or the full Navier-Stokes 
equations. The thin-layer equations do not predict the rise 
in the pressure at the wall because they do not properly 
represent the mechanics of the fully-merged- layer flow. The 
two curves drawn for the full Navier-Stokes equations, with 
and without slip boundary conditions, indicate a rise in 
pressure consistent with the experiments. 
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This clearly shows the superiority of the full Navier-Stokes 
equations over the approximate equations in predicting the 
flow "behaviour in the rarefied regime. 

Finally, in Figure 4*22, we compare the skin-friction 
variation predicted by the thin-layer equations with the 
results of the full Navier-Stokes. equations. The curves are 
drawn for the slip and no-slip cases. Here also, we observe 
the remarkable deviation of the thin-layer solution from the 
solution of the full Navier-Stokes equations at low Reynolds 
numbers. At Reynolds numbers greater than 100, the curves 
of the full Havier-Stokes equations and the thin-layer 
equations do not show any difference. However, as the Reynolds 
number decreases the thin-layer equations predict lower values 
of the skin friction than the full solutions. The difference 
increases as Re decreases. Further, we observe that within 

the frame work of both sets of equations, slip boundary 
conditions increase the skin friction. Also , this increment 
is more for lower values of Re. Finally, we observe that the 
thin- layer equations fail to predict the skin-friction varia- 
tion even qualitatively correctly below Re = 50. 

From the various comparisons made above, it follows 
that the thin-layer equations do not correctly predict the 
phenomenon of merging below Re = 50. For Reynolds numbers 
greater than 50, the thin- layer equations give results with a 
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maximum difference of about 5 percent from the full Ravler— 
Stokes solutions. In brief we can say that low Reynolds 
numbers, the thin-layer equations fail to predict even the 
qualitative behaviour of the pressure and skin-friction 
of the gas at the body. Surprisingly, the thin- layer equations 
give at low Reynolds numbers thicker profiles than the full 
Eavier-Stokes equations. The thin-layer equations over-estimate 
the effect of merging, predict weaker and thicker shocks and 
lower the maximum temperature. It should be noted that all the 
statements made here are based on the cold-wall solutions. 

The same conclusions are valid for the adiabatic -wall case, 
as is evident from the next chapter. Again, it should be men- 
tioned that all the comments made about the thin-layer assum- 
ption are with reference to the thin- layer, equations (2,4-3) 
to (2.47). Stronger forms of thin-layer assumptions have been 
used in literature. For example, Cheng’s (1961, 1963a, 1966) 
thin-shock-layer equations neglect several terms that are 
preserved in the present thin-layer equations. The. conclusions 
mentioned in this section may not be true for the other types 
of thin-layer equations, 

4*5 Effect of Mach Humber; A Case-Study 

In this section, the effect of the Mach number on the 1 
flow is studied at a fixed, typically low Reynolds number, | 




112 


^ e oo = 4-°* Here, Re^ = — — — — , where are the 

Hx> 

free stream density and velocity, r^ the body radius and 
the viscosity at the free stream temperature, - In our usual 
Reynolds number, Re, the viscosity is calculated at the free 
stream stagnation temperature, while in Re^ viscosity is 
calculated at the free stream temperature itself. Calculation 
Re from Re TO depends on f, 1^, and the viscosity temperature 
relationship. In the solutions presented in this section, we 
have used a square-root y-T relationship, and y= 5/3. The body 
temperature is kept equal to the free stream temperature. If 
the free stream temperature is kept constant, increasing 
would mean that the ambient free stream velocity is propor- 
tionally increased, Then, a fixed Re^ implies that with the 
increase of the free stream density is decreased. The 
free stream mean free path increases proportionally and we 
expect the rarefaction effects to become more prominent. 

Numerical integrations are carried out from = 4 
to = 24, at intervals of 2. All the integrations are made 
for the full Navier-Stokes equations with slip boundary 
conditions. 

The tangential-velocity profiles at = 4, 12 and 
20 are shown in figure 4.23. As the Mach number is increased 
the profiles become thicker. 95 percent of the free stream 
value is obtained at r = 1.21 for = 4, at r = 1.36 for 
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^Lo = 12, and at r = 1*5 at = 2C. The slip velocity also 
increases with the Mach nunber. u(l) is 0.15, 0*21, and 0.27 
for 3^ = 4, 12 and 20 respectively. 

Figure 4.24 gives the normal-velocity profiles at 
the Mach numbers of 4, 12 and 20. 95 percent of the free 
stream velocity is obtained at r = 1.45, 1.72 and 2.04 for 
Mbo = 4-, 12 and 20 respectively. Hear the bodies the profiles 
are not much sensitive to the variation of the Mach nunber* 

But the flow is very much affected for r > 1.2. Figures 4*25 
and 4.26 give the temperature and pressure profiles respec- 
tively for the same values of 3^ as above. These figures also 
show a thickening of the profiles with increase of the Mach 
number. The wall pressure p ( 1 ) as well as the ratio fCO/T^ 
increase with the increase of M^. 

In Figure 4.27, we draw the slip velocity, u(l), and 
the temperature, T(l), against the Mach number. We observe 
that, as the Mach number increases, the slip velocity steadily 
increases. However, the rate of increase is small for Mach 
numbers greater than 16. In contrast, the ratio TO)/!^ incre- 
ases rapidly at high Mach numbers. Figure 2.28 shows the 
steady increase of both the skin- friction and heat-transfer 
coefficients with the increase of the Mach number. 
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4.6 Effect of Wall Temperature And Comparison 
With Thin-Two-Layer Theory 

For the purpose of studying the effect of wall tem- 
perature on the flow field at low Reynolds numbers, we take 
a. typical case of Re = 6.8, l^o = 10, Y = 5/3, Pr = .75 and 

solve the full Navi er- Stokes equations (2.15) to (2.21) with 

• ' T 

— w 

slip boundary conditions, for wall temperatures T w = 

0 °° 

ranging from 0.1 to 1.2. At T w =1.2 we almost attain the case 
of an adiabatic wall, with vanishing temperature gradient at 
the body surface. A square-root viscosity-temperature relation- 
ship is used. The present parameters correspond to a fixed 
free stream Reynolds number, Re^, of 40. 

Figure 4.29 shows change in the temperature profiles . 
as the wall temperature is increased. .The profiles are drawn 
for Tyy = 0.2 to 1.2 at intervals of 0.2. We observe that, as 
the wall temperature is increased, the extent of disturbance 
become larger* Eor T w = 0.2, the temperature attains a value 
of T = 0.1 at r * 1 . 67 , while at T w = 1.2 this value is attain- 
ed at r -1.83. With the increase of wall temperature, the 
temperature gradient at the body decreases continuously. 

At T w = 1.2, this gradient almost vanishes, resulting in an .. 
adiabatic-wall profile. Figure 4.30 gives the normal-velocity 
profiles at three different wall temperatures, namely., . . . 

T w = 0.2, 0.6 and 1.0. As in the case of temperature profiles, 
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the velocity profiles also show an increase in the extent 
of disturbance. The variation of the wall temperature does 
not have much effect ‘on the normal velocity upto a certain 
distance from the body surface. However, near the body the 
normal velocity itself is very small. 

In Figure 4.31, we present the variation of the slip 
velocity and the density at the wall with the wall tempera- 
ture. The slip velocity steadily increases with the increase 
of wall temperature* As the wall temperature T w is raised 
from 0.1 to 1.2, the slip velocity is almost doubled. In the 
same range, the density p~ of the gas near the body decreases 
from about 13.0 to 3.5. The density decreases faster at low 
wall- temperatures than at high wall-temperatures. With the 
decrease of density, we expect the slip effects to become 
prominent. This is reflected in the increase of slip velocity. 
However, the temperature jump decreases with the increase of 
wall temperature owing to the decrease in the temperature 
gradient at the body. The gradient of the tangential velocity 
also decreases, because of the increase in the extent of dis- 
turbance as the wall temperature increases. But this decrease 
in the tangential-velocity gradient is not large enou^i to 
outweigh the effect of the increase of the mean free path 
near the body. The net result is an increase in the slip 


velocity 



116 


In figures 4.32 and. 4.33, the gas temperature at the 
wall and the heat-transfer rate to the "body are drawn against 
the wall temperature. Corresponding curves are also drawn 
for the thin-layer theory of Mu (1967). Mu’s results are ' 
based on the formulation of Cheng, discussed in detail in 
Section 1.4* Mu treats the effect of slip velocity and tem- 
perature jump as a perturbation on the no-slip solution, his 
perturbation parameter being . The no-slip solutions 

of this scheme are exactly the same as the basic solutions 
of Cheng (1961), developed for the incipient-merged-layer 

regime. The results of these theories are given in terms 

2 

of a rarefaction parameter, K , which is related to Re as 
follows : 
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where Re is the Reynolds number given in the equations 

(2.14), and the subscript * refers to an appropriate reference 

temperature T*. To compare our results with those of Mu, It 

2 

is necessary to obtain the K corresponding to our solutions. 

The reference temperature, suggested by Cheng and Chang (1963') 

T + T 

as being representative of the flow, is T* = , where 

T s is the temperature behind the shock wave as determined by 
the thin— shock— layer theory , which takes into account the 
shock-slip conditions. Depending on the rarefaction of the 
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flow 5 this T g is less than the temperature predicted by the 
Rankine-Hugoniot shock conditions. The more rarefied in the 
flow, the more T g deviates from the Rankine-Hugoniot value. 
Since the concept of ’behind-the-shock-wave * does not enter 
naturally into our present numerical solutions, where the 
shock transition comes as a part of the solution, we take 
the maximum temperature T_ _ in the flow field instead of I 

Illcb'X. £3 

to determine the rarefaction parameter ■. Thus, different 
2 

K s are obtained for different wall temperatures. It is 

observed that within this variation of Liu’s results do 

not show any appreciable change. The results of Liu presented 

2 

in figures 4.32 and 4.33 are for an average value of K = 1.35. 

In figure 4.32, we present the variation of temperature 
at the wall with slip boundary conditions with prescribed 
wall temperature, The straight line in the figure corresponds 
to the case of neglecting the slip boundary conditions. The 
temperature jump predicted by Liu’s theory is always less 
than that of the present solutions. The temperature jump is 
given by the height of the corresponding curve above the 
straight line representing the no-slip case. In most of the 
region, Liu’s prediction is about half of the prediction of 
the full Navier-Stokes equations. The figure shows that for 
the full Mavier-Stokes equations the temperature jump vani- 
shes at ? * 1 . 2 . Since the temperature jump vanishes when 
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T ' at the body surface is zero, this temperature is the adia- 
batic recovery temperature. In the case of Liu's theory, the 
temperature jump always vanishes at T = 1, irrespective of 
the rarefaction parameter K . Experiments (for example, Drake 
and Backer, 1952, Hickman and Giedt, 1963) show that in rarefied 
flow the recovery temperature is higher than the free stream 
stagnation temperature. Thus, Liu’s theory fails to predict 
this experimentally observed phenomenon. 

Figure 4.33 gives the variation of the heat-transfer 
rate to the body with wall temperature. Our results shown 
are from the solutions, of the full Navier-Stokes equations 
with slip boundary conditions. For Liu’s theory, the slip as 
well as no-slip curves are drawn. Liu's slip-correction to the 
heat-transfer rate is calculated from the following formula 
(from equation (5.1) of Liu's (1967) paper). 
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C TJ is the no-slip heat-transfer rate given by 
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where k = 


2 Pr 


{ / -K + 1 - 1 ) 

/ g/ 


(4.4) 


^ 1 

\ ( 3*5 k) is the incomplete gamma function of order 1/5 and 

is related to the integral 


a V3 


I ( a 1 / 5 ) = J P(i a) = f exp(-t 5 ) dt 
^ J 0 


(4.5) 


This integral has been directly calculated and tabulated by 
Abramowitz (1951). Similar slip-correction to heat transfer 
is also obtained by Cheng and Ghang ( 1965 ). 


Figure 4.55 shows that the heat- transfer rate pre- 
dicted by Liu's theory is always less than that predicted 
by the full Eavier-Stokes equations* The slip-correction, 
given by Liu is very small and always reduces the no -slip 
heat transfer rate. Even though we do not calculate the no- 
slip results for the full Navi er- Stokes equations .for the 
entire range of wall temperatures, the single solution at 
T w = 0.029 (please see Figure 4*15) shows that the slip- 
correction to the heat transfer rate is not small at Re = 6.8 
(which is the Reynolds number for the present case). Moreover, 
the effect of slip conditions is not to decrease the heat- 
transfer rate as predicted by Liu's theory, but to imcroase- 
it.- This disertpancy may be because of the fact that Liu's 
theory underestimates the temperature jump at the body, as is 
already seen in Figure 4 . 52 . 
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4.7 Comparison With Other Theories 

In the previous section, we have compared our results 
with those of Liu (1967), with reference to the effect of 
wall temperature on the flow at one particular Reynolds 
number. In the present section, we make further comparisons 
with the results of the thin- two-layer theories and also 
with some results of Monte Carlo simulation for highly 
rarefied flows. 

In figure 4.34, we draw the slip velocity and tempe- 
rature jump versus K 2 , as obtained from the full Ravier- 
Stokes equations for the values of the parameters, y = 1.33,. 

Moo = 10, Pr = 0.71 and 3^ =0.1, and compare them with corr- 
esponding results from the work of Liu (1967). The rarefac- 

2 

tion parameter K is determined from Re in the manner already 

described in the previous section, These curves show that 

the temperature jump and slip velocity predicted by Liu’s 

theory are about 50 percent less than those predicted, by the 

full Navier-Stokes equations in almost the entire range of 
2 

II presented in the figure. This underestimation of the slip 
quantities by Liu is due to the a priori assumption that the 
slip effects are small and treating them as perturbations to 
the no-slip solution. Purther, these effects are calculated 
within the framework of the thin-layer assumption, which 
itself is not valid in highly rarefied flows. It may be noted 
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that according to the approximations inherent in liu's 

theory, the mean-free-path x w , appearing in the equations 

(2.8) and (2.9) for the slip velocity and temperature jump, 

is calculated at the body temperature T , rather than at the 

w ■ 

temperature of the gas near the body (see equation 2.18 of 
liu‘s paper (1967)). This because of the particular pertur- 
bation scheme adopted by liu. 

In Figure 4.35, we compare the temperature profile 
calculated from the full Navier-Stokes equations, for the 
parameters Re 0 = 57.4 , y= 1.4, Pr = 0.75 and - 18.8, with 
a corresponding experimental profile measured by ihouse and 
Bogdonoff (1969) and the profile calculated from Cheng's 
theory. The experimental curve is given as a shaded band in 
accordance with the estimated experimental accuracy. The curve 
given for Cheng's theory is also reproduced from the paper 
by Ahouse and Bogdonoff. This figure clearly shows the failure 
of Cheng's two-thin-layer-model theory in predicting the 
detailed flow structure, even though it predicts the heat- 
transfer rates to the body in a good agreement with experiments. 
Cheng's profile shows a larger maximum temperature and a thi- 
nner shock wave than the experimental profile. The total 
extent of disturbance is nearly half of what the experiment 
shows. €tn the other hand, our profile calculated from the full 
Navier-St okes equations agrees reasonably well with the experi- 
mental profile. Our profile also predicts a slightly higher 
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maximum temperature than the experimental profile, but this 
not appreciable in comparison with the large difference 
shewn by Cheng’s theory* 

Comparing the profile obtained from Cheng’s model 
with their experimental profile, Ahouse and Bogdonoff have 
conjectured that a Navier-Stokes approach is inadequate to 
analyse the phenomenon of merging. However, the present resu- 
lts show that a formulation in terms of the full Ravi er- Stokes 
equations gives a good agreement with experiments. 

In Bigure 4*36, we compare our normal -velocity and 
temperature profiles with those obtained from Cheng’s two- 
thin-layer model at He = 70 and 26.6. The profiles for Cheng’s 
theory are calculated by Mu and S 0 game (1969) f or y = 1.4 and 
a very cold wall. Our profiles are obtained from solving 
the f ull Navier-Stokes equations with y = 1.4, = 20 and 

5^ = 0.05 and using a linear viscosity-temperature relation- 
ship, which is implicit in Cheng’s theory. One notes that the 
profiles of normal velocity and temperature given by Cheng’s 
theory are discontinuous at the ’shock interface’. These 
discontinuities are unavoidable . in Cheng’s two-layer formula- 
tion which matches the derivative of v at the interface. At 
Re = 70, it is observed that our profiles and Cheng’s, profiles 
agree well throughout the layer of disturbance.. However, .at 
Re = 26.6, Cheng’s profiles differ very much from our profiles 
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in the shock transition-zone, while the agreement is good 
from the body surface. to the shock interface. Cheng's solu- 
tions underestimate the shock thickness* 

A comparison with Cheng's theory at more rarefied, 
conditions is presented in figure 4.37. Here, . Cheng's profile 
(as given by liu and Sogame, 1969) is for K 2 = 1 , In our 
computations we have used a square-root viscosity-temperature 
relation. Since Cheng's theory uses a linear viscosity law, 
the transformation of K 2 to Re should take into account a. 
ref erence-temperature concept as given in the Equation (4*1). 
We calculate Re using the ref erence-temp erasure concept 
(which gives Re = 11.5j> as well as neglecting it (which gives 
Re = 7.0), Agreement between Cheng's profile and- our .profile 
for Re = 7.0 is good near the |body surface upto r -1.2. How- 
ever the shock-transition zone predicted by Cheng’s theory 
is much thinner compared to our profile. This difference in 
the extent of disturbance is lessened. if we compare Cheng’s 
profile with our profile at Re = 11.5. But, the agreement of 
the flow near the body is lost. Further, the shock-transition 
zone given by our profile at Re = 11.5 also is thicker than 
Cheng’s prediction. Thus the comparisons in Figures 4.36 and 

4.37 indicate that Cheng's thin- two -layer theory predicts 
thinner layers of disturbance than the full Navier— Stokes 
prediction and that the differences increase as the Reynolds 
number decreases.. 
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In Figure 4.38, we present a comparison of the full 
Navier-Stokes solution at Re^ = 152 with the results obtain- 
ed by Vogenitz and Takata (1970) using Monte Carlo direct 
simulation technique, Yogenitz and Takata calculate the flow 
in the nose region of a spherically blunted slender cone 
(8° half -angle). Collisions of the molecules with the body 
and with each other are accounted for, The reflection from 

the body is assumed to be diffuse with perfect thermal acco- 
mmodation (cr = a = l). In our integrations of the full 
Navier-Stokes equations with slip boundary conditions, we 
use a linear viscosity-temperature law and unit accommodation 
coefficients. Figure 4.38, where the normal-velocity profiles 
are compared, shows that our solution agrees well with the 
Monte Carlo profile in the entire range disturbance, 

4.8 Comparison with Experiments 

In this section, we present detailed comparison of 
our results obtained from solving the full Navier-Stokes 
equations (2,15) to (2.21) with the available experimental 
results. In Section 1.6, we have already reviewed briefly the 
experimental results for the blunt-body problem under hyper- 
sonic rarefied conditions. Here, we shall compare our theore- 
tical predictions with the experimental data for the stagna- 
tion-point heat -transfer rates, the temperature profiles and 
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the density profiles* 

Comparison With Stagnation-Point Heat-lransf er Rates { 

The heat-transfer results are given in Figures 4.39 
and 4*40. In Figure 4.59 we compare our results with the 
experimental results of Ferri, Zakkay and Ting {1962% The 
heat -tr ansf er rate is presented in terms of the parameter 

(v 
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where H is the total enthalpy. This Gg is presented 
against the Reynolds number Re^, defined by 
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The experimental results are available in terms of these 

parameters, and we have transformed our results into these 

new parameters. In the present case, the Reynolds number 
p V r -q 

Re = — — — — is related to Re^ by 

v 

0 °° 

Re-™ = 5.207 Re 

i* 

The experiments of Ferri, Zakkay and .Ting (1962), shown 
in the figure, were carried out at a free stream stagnation 
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temperature of 2100°R and the corresponding value ofy sugges- 
ted was y= 1.326. Figure 4*38 shows a very good a.greement 
of our results with the experimental data in the entire 
range of the Reynolds number. 

In Figure 4.40, we make comparisons with some other 
experimental heat— transf er data, some of which extend to 
even more rarefied conditions than those in Figure 4.39* 

The experimental data are taken from Ferri, Zakkay and Ting 
(1961), Wittliff and Wilson (1962), Vidal and Witfcliff (1963)., 
and Bovlan (1971). As indicated in the figure, the conditions 
at which these various experiments are performed are much 
different from each other. The theoretical curve predicted by 
Cheng’s thin-layer theory and the boundary-layer theory are. 
shown in the figure. Also shown is a curve based on a theory 
for vo rticity- interact ion regime {Cheng, 1961), using the 
usual Rankine-Hugoniot shock relations. The vorti city-inter- 
act ion curve deviates upwards from the experimental data 
for K 2 less than about 3. The agreement of our results as 
well as the prediction of Cheng’s thin-layer theory with 
the experimental data is reasonably good. In view of the 
differences in various essential parameters, it is not 
possible to make any conclusion regarding which theory 
shows a better agreement. The heat-transfer results show 
an appreciable dependence on the ratio of the specific 
heats Y, or in other terms, on the density ratio across a 
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normal shock, This point is clearly shown in Figure 4*1 5, 
where the heat transfer rates at Y = 1,667 are compared with 
those at Y = 1,33. This importance ofy is also observed by 

Potter (1963) and Boylan (1971). Potter suggests a heuristic 
correction for the heat-transfer rates when results are to 
be compared corresponding to different values of Y. However, 
we do not apply any such correction, but present the experi- 
mental and theoretical results as they are for various ys. 

It may be observed from Figure 4.40 that at high 
values of the parameter E^, various theoretical curves tend 
to the corresponding boundary-layer limit shown in the 
figure. It is of interest to note that our results cross the 
value unity at about = 0.2. A unit value of Cg corresponds 
to the value given by the free-molecul e-flow theory for unit 
accommodation and reflection coefficients. Such an overshoot 
of the heat transfer over the free-molecule limit is also 
predicted by the Monte Carlo simulation studies of Bird 
(1966). His results, at a molecular speed ratio of 10, show 
that for a highly cooled body, the heat transfer rises above 
the free molecule limit at a free stream Khudsen number of 
about 1.0. At Khudsen numbers about 30 the heat -transfer, 
rate again decreases smoothly to the free-molecule value. 

The Khudsen number Kn = X^r-g at which our heat-transfer 
results cross the free-molecule value is also of the order of 
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unity. However, further decreasing the Reynolds number, our 
heat-transfer results go on increasing and do not show any 
tendency to return to the free-molecule limit. This indicates 

that for values of Eh greater than unity, the full Favier- 
Stckes equations fail to give correct results. It may he 
noted that Bird*s predicted overshoot of the heat transfer 
above the free-molecule value is neither verified nor disp* 
roved by the available experimental heat-transf er results* 

T empera.turc-p.rof ile comparisons ; - 

In the Figures 4.41 to 4.44 we compare some solutions 
of the full Havier-Stokes equations with slip boundary, condr 
itions with the experimental flow-field studies made by Ahous 
and Bogdonoff (1969). In the Section 4.7, we have already 
compared our profile at Re^in = 5060 with the corresponding 
experimental profile. In the present section we compare four 
more of our temperature profiles with experiments, for the 
conditions given in Table 4*1. Ihe last row of the table 
corresponds to the comparison made in Section 4.7. 
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TABLE 4.1 - Test-Conditions for the Experimental Temperature' 
Profiles, With Which Comparisons Are Made 


Re-oo /in 

K> 

V(°K> 

Re 

s 

Re 'Re 

using Suther- using linear 

land law viscosity law 

40780 

18.9 

1600 

461 .0 

316.2 

212.3 

15400 

26.7 

1800 

91.2 

51.48 

40.22 

7550 

24.6 

1800 

51.4 

31.67 

23.20 

7420 

20.8 

1600 

68.8 

39.81 

31.79 

5060 

18.8 

1600 

57.4 

39.81 

26.47 


The first four columns are taken from the test- 
conditions given by Ahouse and Bogdonoff. Ro g , the Reynolds 
number behind a normal shock wave, is obtained from comple- 
mentary measurements made upstream and downstream of a normal 
shock. The last two columns are obtained by transforming the 
free stream Reynolds number, Re^, to the Reynolds number Re, 
by Sutherland and linear viscosity-temperature laws respectiv- 
ely. 

Ahouse and Bogdonoff have measured the rotational 
temperature profiles along the stagnation stream line of a 
three-quarter inch diameter, water-cooled, stainless steel 
hemisphere-cylinder, in a nitrogen stream, using electron 
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beam fluorescence technique. In the measurements, a wall to 
free stream stagnation temperature ratio of about 0, 19 . is 
maintained. We compare our thermodynamic temperature profiles 
with the experimental rotational -temperature profiles* This 
is permissible because the rotational degrees of freedom for 
nitrogen are excited at temperatures as low as 2.9°K 

(Zel'dovich and Raizer , 1966), and it is always possible to 
assume that the rotational temperatures are equal to the 
translational temperatures, 

im important parameter that one should be cautious 
about is the Reynolds number. Basically, in all the present 
numerical integrations we use the Reynolds number Re defined 
Voo r B / u 0 oo> which comes as a result of our nondimen- 
sionalisation. The experimental conditions are given in terms 
of the free stream Reynolds number Re^. The change of Re^ to 
Re requires the knowledge of proper viscosity-temperature 
relationship. In general, in all our numerical integrations 
wc have used the square-root viscosity- temperature relation- 
ship. However, if we change Re TO to Re using wa fT , it is 
found that the corresponding Re is too big and the solution 
gives a profile much thinner than the experimental profile* 
This is because a realistic viscosity-temperature relation- 
ship is not used. This can be seen from the fact that the 

measured Re (the Reynolds number behind a normal shock wave), 
s 
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as shown in the column 4 of Table 4.1, differs very much 
from, the calculated Re g , using square-root viscosity law. 

It is proper to calculate Re from a more realistic u-T rela- 
tionship, like the Sutherland law. Even a linear viscosity- 
temperature relationship is found to be good. Columns 5 and 6 
of Table 4*1 give the value of Re obtained using the Suther- 
land and linear viscosity laws. Our numerical integrations 
show that, once Re is calculated, the solutions obtained are 
not very much sensitive to the actual v-T relation used in 
the iterative procedure, at least in the moderately high. 
Reynolds number cases studied in this section. This suggests 
that Re = Pqo ^ r-g/w^, or any other Reynolds number for 
which the viscosity is measured at a temperature comparable 
to is a far better correlation parameter for the rarefied 

blunt -body problem, than the free stream Reynolds number Re^ 
or a comparable Reynolds number. 

In the comparisons presented in figures 4.41 to 4.44, 
the experimental profiles are drawn as thick bands in accor- 
dance with the estimated experimental accuracy reported by 
Ahouse and Bogdonoff. figure 4.41 shows, that, in the high 
Reynolds number case of Re^/in. = 40780, our full Navi er- 
Stokes results compare reasonably well with the experiments. 
Our profile, however, shows a slightly sharper shock than, 
the experimental profile-, figures 4.4-2, 4.43 and 4.44 show 



that the full Ravier-Stokes equations show a good agreement 
’with the experimental results, except near the peak tempera- 
ture region. In all the cases the peak temperatures predicted 
by the full Ravier-Stokes equations are slightly higher than 
the experimental values* In all the figures except in 
Figure 4*44, we present the solutions for which Re is calcu- 
lated using linear viscosity-temperature relationship. In 
Figure 4*44, we show, in addition, the profile .obtained by 
calculating Re by the Sutherland viscosity law. It can be 
observed that the difference in the two solutions is small. 

It may be noted in all these figures, that, in a 
region near the body surface the experimental profiles are 
not drawn. This is because of the difficulty in measuring 
the temperatures in this region due to significant three-dimen- 
sional effects. This precludes any experimental conclusions 
with regard to the temperature jump* 

Density-profile comparisons 

Russel (1968) has measured the density profiles ahead 
of a sphere in rarefied streams of argon and nitrogen, using 
an electron-beam, X-ray technique. In Figures 4.45 to 4.49, 
we compare our profiles, obtained from the full Navier— Stokes 

i 

equations using slip boundary conditions, with Russel’s 
experimental profiles. Ihe conditions, at which these 
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comparisons are made, are presented in Table 4,2. 


Table 4.2. - Test -Conditio ns for the Gold-Wall, Bxperi- 


are Made. 


M m 

CO 

T 

w 

V ( ° K > 

Re „ 

s 

Re (using linear 
viscosity law) 




Argon 



100 

3.83 

.26 

300 

19;8 

16.97 

200 

3.83 

.26 

300 

39,5 

33.94- 

500 

3.80 

*26 

300 

91.5 

85.97 




Nitrogen 


30 

4.20 

.26 

300 

7.7 

6.62 

500 

4.11 

.26 

300 

134.0 

114.2 


In our numerical integrations we use the values of 
Re as given in the last column of Table 4.2, which are obtain- 
ed from Re^o using linear viscosity-temperature law. The com- 
parisons include very low and higb Reynolds number cases, 
figures 4.45 to 4.49 show a good agreement of the present 
solutions with the experimental profiles. Our results show 
a slightly thinner extent of density— disturbance than the 
experimental profiles. However, the agreement is better in 
low-Reynolds-number cases. 


V, u 

T, P 



FIG. 4.1 _COLD-WALL FLOW PROFILES; 
Re = 296 , V = 5/3 



FIG.4.2_C0LD-WAI_L FLOW PROFILES 
Re = 34 , V = 5/3 
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FIG. 4.3- COLD -WALL FLOW PROFILES , Re = 3 .4 , V = 5/3 



NUMBERS, COLD WALL , ? = 5/3 
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FULL NAVIER- STOKES EQUATIONS 
SLIP CONDITIONS 

Moo = 10 

T =5/3 
Pr= 3/4 


FIG. 4.5_ TANGENTIAL- VELOCITY PROFILES 
AT VARIOUS REYNOLDS NUMBERS 
COLD WALL , 1 =5/3 
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FIG 4 8- TEMPERATURE PROFILES AT VARIOUS 
NUMBERS; COLD WALL , V =1.33 
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FIG. 4. 9_ PRESSURE PROFILES AT VARIOUS REYNOLDS 
NUMBERS ; COLD WALL , V =1 33 
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FIG. 4.10. TANGENTIAL- VELOCITY PROFILES AT VARIOUS 
REYNOLDS NUMBERS ; COLD WALL , V =1.33- 




FIG. 4.12 -COMPARISON WiTt- 
AT Re =3.4, 1 = 5/3 
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FIG.4. 15.EFFECT OF SLIP BOUNDARY CC 
ON HEAT - TRANSFER COEFFICIENT 
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.MAXIMUM TEMPERATURE IN THE FLOW FIELD VERSUS REYNOLDS 
NUMBER . FULL NAVIER - STOKES AND THIN • LAYER EQUATIONS. 
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FIG. 4.20- EXTENT OF DISTURBANCE VERSUS REYNOLDS 
NUMBER, FULL NAVIER- STOKES AND THIN - 
LAYER SOLUTIONS 
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FIG. 4.26 












9k 



FULL NAVIER 






VToco 

FIG. 4.31 -SLIP VELOCITY AND DENSITY AT 
THE BODY SURFACE Vs WALL 
RATURE 
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FIG. 4.32-TEMPERATURE JUMP Vs WALL TEMPERATURE 
COMPARISON WITH LIU'S (1967) THEORY 




FIG. 4.33. STAGNATION- POINT HEAT- TRANSFER RATE 
Vs WALL TEMPERATURE, COMPARISON W! 
LIU (1967) 
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CHENG'S THEORY 
Re s =60> i =1.4 

Twall/T 000 = 0,2 . 

PRESENT SOLUTION 
Re s =57.4, 1=1.4 
T V1IAI ,/T =0.19 

'WALL ooo 

EXPERIMENT 

/ A HOUSE & BOGDONOFFCI969) 
X NITROGEN 
/>\ Rec = 57.4 , 


FIG. 4.35-COMPARISON OF T Pi^™ R T E HE oRY F ' LE 
EXPERIMENT AND CHENG S THEQR 


WITH 


PRESENT 



PROFILE COMPARISON WITH CHENG'S 

TWO -LAYER MODEL (FROM LIU & SOGAME,1969) 

AT Re = 70 & 26.6 







FIG. 4.38-COMPARISON OF NORMAL - VEL 
MONTF CARLO CALCULA TIONS 
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FIG 4.42 -COMPARISON OF TEMPERATURE PROFILE 
WITH EXPERIMENT. NITROGEN , Re^ 5775 
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FIG. 4.4 5 -DENSITY - PROFILE COMPARISON WITH 

EXPERIMENT , ARGON , Re s = 19.8 > COLD WALL 
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CHAPTER 5 

RESULTS FOR THE AD LAB .AT T C-W h T.T. CASE 

5,1 I ntroduction 

Solutions of the governing equations with adiabatic- 
wall condition are discussed in this chapter. The full 
Navier-Stokes equations (2,15) to (2,21) and the t hin -layer 
equations (2.43) to (2,47) are solved with the adiabatic- 
wall condition, namely, !P (1) = 0. Solutions are. sought 
with slip as well as no-slip boundary conditions. Since the 
temperature gradient .is zero at the body surface, the tem- 
perature jump given by (2,32) always vanishes. The tempera- 
ture of the gas near the body is equal to the temperature 
of the body itself. But, this temperature is not known a 
priori and is obtained as a part of the solution. 

The main numerical integrations performed are for 
the case of a monatomic gas stream. The full Navier-Stokes 
equations are solved for the following parameters; 

0 ) Y - 5/3, K, = 10, Pr = 3/4, 

(2) 2 < Re < 300, 

(3) v = y5 

(4) o =1.0 


9 
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Several solutions of the thin-layer equations also are 
obtained in this range aprapos determining the validity of 
these equations. In order to make comparisons with experi- 
mental density profiles of Russel (1968), and the pressure 
data of Sherman (1953) and Potter and Bailey (1963) several 
integrations are performed at various Reynolds numbers and 
Mach numbers for monatomic and diatomic gas streams. 

In the next section we present some general features 
of the adiabatic wall profiles. The effect of slip velocity 
is investigated in Section 5.3. Section 5.4 contains a dis- 
cussion about the validity of the thin-layer equations with 
reference to the adiabatic-wall solutions. In Section 5*5 
we compare our results with the predictions of the second- 
order boundary-lay or theory. In Section 5.6 we make a detail- 
ed comparison of our results with the experimental data 
for density profiles and impact pressures. 

5 • 2 de n oral Plow -Prof il es 

The; general structure of the adiabatic-wall solutions 
is presented in this section through figures 5.1 to 5.8. 

All the results discussed in this section are obtained 
from the full Navier-Stokes equations with the slip boundary 
condition. The integrations are performed at = 10, 
y = 5/3 , Pr = 3/4 and o~ = I* 0 * 



In Figure 5.1, we draw the velocity, temperature ana 
density profiles at Re = 170. The total extent of disturbance 
caused by the body is only about 23*/. of the body radius. There 
is a quick transition in the values of the normal velocity, 
temperature and density in the range 1.19 < r < 1.22, which 
represents the shock wave. Between the body and the inner 
edge of this region, the variation of temperature, density and 
hence of pressure is very small. The slip velocity at. the body 
and at r = 1.06 attains 80/. of its free-stream value. At this 
distance the normal velocity attains only 5/. its free-stream 
value. These factors indicate that the usual boundary-layer 
assumptions can be considered to be valid at this Reynolds 
number. 

In Figure 5.2, the profiles obtained at Re = 34 are 
presented. The extent of disturbance increases by about 75/. 
compared to the case of Re * 170. Also, a steep shock-transi- 
tion zone does not exist in this case. The shock wave thickens . 
appreciably and merges with the viscous layer near the surface. 
This merging results in an appreciable gradient of pressure, 
and density near the body surface. The boundary-layer assum- 
ptions are not valid in this case. 

The profiles drawn in Figure 5.3 are at Re = 4.25 and 
correspond to veiy rarefied conditions. Wo trace of a 
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shock-transition zone exists. The extent of disturbance is 
very large and .is even greater than the body radius. The 
slip velocity at the body is as much as half of the free- 
stream value of the tangential velocity. Even the noimal 
velocity, v, changes rapidly near the body* The density and 
hence the pressure have very steep gradients near the body 
surface. At moderately high Reynolds number flows past an 
insulated body, it is usual to assume a constant density in 
the region between the shock wave and the body. Eor example, 
Probet cin and Kemp (i960) use this constant-density assump- 
tion in their viscous-and incipient-merged- layer studies and 
Oberai (1964) uses it in his third-order boundary- layer 
calculations. However, such an assumption is not valid in the 
low-Reyno Ids ^number cases as is evident from Figure 5.3. 

In Figures 5.4 and 5.5, we draw the profiles of the 
tangential velocity and the normal velocity respectively . at 
five different Reynolds numbers, ranging from the very low 
to high Reynolds number cases. The tangential velocity pro- 
files show an increase in the slip velocity with decreasing 
Reynolds numbers. It may be noted that, at a given Re, the 
tangential velocity attains its free-stream value much nearer 
to the body than .the normal velocity. Thus, the tangential- 
velocity profiles do not give a correct indication of the 
extent of disturbance caused by the body. The normal -velocity 
profiles clearly indicate the thickening of the shock-wave -like 
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region. However, upto a certain distance from the body, the 
change in uhe normal velocity with Reynolds number remains 
r olat iv oly small . 

In Figure 5.6, we present the density profiles. The 
ratio o.t the density at the wall to the free— stream density 
doerea;;;',!.; \.ith the decrease of Reynolds number. Also, the 
structure of the profiles near the body changes appreciably. 

At low Reynolds numbers, there is a steep variation of den- 
sity near the wall. At high Reynolds numbers this variation 
isnegligi’blo. 

The temperature profiles drawn in Figure 5.7 clearly 
show the effect of merging on the flow near the body. The 
profiles are drawn at Re = 3.4, 8.5, 17.0, 85.0 and 255.0. 

The wall temperature shows a continuous increase as the 
Reynolds number decreases. At low Reynolds numbers this 
increase is faster than at high Reynolds numbers. The incr- 
ease of the adiabatic recovery temperature at low Reynolds num- 
bers is mi experimentally observed phenomenon (Drake and 
Baker, 1952, Hickman and Giedt, 1963) cand is related to the 
rise in the pressure at the wall. 

Figure 5.8 shows the variation of the stagnation-point 
pressure with Reynolds number. In addition to the adiabatic- 
wall case the curves corresponding to the cold-wall cases 
also arc drawn. As the Reynolds number is decreased the 
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pressure decreases at first and later steeply increases, 
this is in conformity with the general pattern of the experi- 
mental results (Figure 1*3)* Direct comparison of our results 
with the experimental data of Sherman (1953) and Potter and 
Bailey (1963) is made in a subsequent section and the agree- 
ment is found to be good. The rise in the impact pressures 
and the adiabatic recovery temperatures at low Reynolds num- 
bers is due to the removal of the shielding effect of the 
shock wave which exists at moderately high Reynolds number 
cases. Also, from Figure 5«S we observe that as the wall 
temperature increases, the wall pressure decreases. This is 
in agreement with the results of Potter and Bailey (1963). 

5 • 3 E ffect of Slip Boundary Conditions 

In the adiabatic-wall case the temperature jump always 
vanishes owing to the condition that the temperature gradient 
at the wall is zero. However, the slip velocity exists and 
is larger in magnitude compared to a cold -wall case. This is 
so because near an adiabatic wall the gas is much more rare- 
fied than near a cold wall. 

ft r e find that at Re = 170 there is virtually no diff- 
erence between the slip and no-slip profiles, except in a 
small region near the body. The slip velocity is 0.093. The 
no-slip tangential-velocity profile catches up with the 
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slip piofile at about r — 1.01 after which, no appreciable 
difference exists between the slip and no-slip cases. She 
slip profiles at Re = 170 are shown in figure 5.1. 

In Figure 5.9, the slip and no-slip profiles at Re=17 
are compared. Here we find that slip affects the entire flow 
field from the body surface to the freestream. The slip 
velocity, u(l), is .279. Slip has the effect of decreasing 
the pressure and density levels throughout the flow field. 
However, the temperature (not shown in the figure) for the 
slip case is always slightly higher than the no-slip case. 

The difference is maximum at the body surface, where the 
slip solution gives T = 1.119 and the no-slip solution gives 
T = 1 ,. 090 . 

In Figure 5.10, we present the slip and no-slip 
solutions at Re = 5.1. We find that the effect of slip 
throughout the flow field is further enhanced. The slip velo- 
city is .475. Introducing the slip condition reduces the 
density at the body surface from 5.801 to 2.804. The temper- 
ature T of the wall increases from 1.514 to 1.402. 

In Figure 5.11, we show the variation of the slip 
velocity with Reynolds number. At high Reynolds number the 
slip velocity is very small and tends to zero. At Re = 250, 
the slip velocity is about 8 percent of the free-stream value 
of the tangential velocity. As the Reynolds number is 
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decreased, the slip velocity increases rapidly, and at 

Re = 2, it is more than 50 percent of the free-stream value. 

In the same figure, we show a corresponding cold— wall slip- 
velocity curve. The slip velocity on the adiabatic wall is 
much larger than on the cold wall. 

In iigure 5.12, the extent of disturbance for the 
slip and no-slip cases is drawn against the Reynolds number. 

As in the cold-wall case, the extent of disturb an ce is arbi- 
trarily defined as the nondimensional distance from the body 
to that point where T = 0.1, The figure shows that the total 
extent of disturbance is 'slightly decreased by introducing 
the slip boundary condition. In the cold-wall case, as obse- 
rved in Figures 4.12 and 4.14, the trend is in an opposite 
direction, where it is found that the slip boundary conditions 
increase the total extent of disturbance. This can be expla- 
ined from the fact that in a cold-wall case the temperature 
jump reduces the cooling effect of the body. Cooling of the 
body results in a shrinkage of the extent of disturbance. 
Hcnoc, if the temperature jump is taken into account, this 
shrinkage is reduced. In the insulated-wall case, the tempe- 
rature jump is zero and Figure 5.12 shows that the effect 
of slip velocity alone is to decrease the extent of distur- 
bance. In the cold-wall cases studied in Chapter 4, the 
effect of temperature jump has dominated the opposite effect 
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of slip velocity. 

figure 5.13 gives the variation of skin friction for 
the adiabatic— wall case with slip and no— slip boundary 
conditions. ho observe that the slip velocity decreases the 
skin friction, for Reynolds numbers above 50, the dif ference 
is small. But below this value the difference becomes larger 
and increases with decrease of Reynolds number. Again, the 
effect of slip on skin friction for an adiabatic wall is in 
an opposite direction to that observed for a cold body 
(of. figure 4.16). The temperature jump in the cold-wall case 
increases the coefficient of viscosity jT, which contributes 
to the increase of the skin friction. The effect of the slip 
velocity alone is to decrease the skin friction, as the tan- 
gential-velocity gradient at the body decreases with slip. 

finally, the adiabatic wall temperatures predicted 
by the slip and no-slip solutions are not the same. This is 
shown in figure 5.14, where the wall temperature is drawn 
against the Reynolds number. Solutions taking slip velocity 
into account predict a higher wall temperature than the no- 
flip solutions. However, the temperature jump in both cases 
is zero, and the wall and the gas near the wall have the 
s cm e t emp e r a t ur o , 

Briefly, we have drawn the following important . 
conclusions regarding the slip effects in the stagnation 
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region of an adiabatic blunt body: 

(1 ) Slip exercises an important role in determining 
the flow field* At low Reynolds numbers slip affects the 
entire flow structure from the body to the free stream. 

(2) Compared to a cold-wall case, the effect of slip 
conditions on the flow is more in an adiabatic-wall case. 

(3) Slip decreases the skin friction and slightly 
lessons the extent of disturbance in the gas. 

(4) She nature of slip effects depends on whether the 
body is cold or hot. The slip velocity and the temperature 
jump have opposite effects on certain quantities like skin 
friction, and which dominates the other depends on the tem- 
perature of the body, 

5.4 Limit of Validity of the Thin-layer Assumption 

In solving the thin-layer equations (2.43) to (2.47) 
for the adiabatic-wall case, either with slip velocity or 
neglecting it, the Becker's integral (2.42) is used. This 
saves a considerable amount of the computer time. The use of 
Becker's integral in the present case is made possible bec- 
ause we take the Prandtl number equal to 0.75. A value of 
Pr =0,667 is consistent with the monatomic gas studied here. 
However, it is verified in some individual cases of a cold 
body as well as an insulated body, .th at this difference in 
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the Prandtl number does not appreciably alter the solutions 
obtained* Uoing a smaller Prandtl number results only in a 
small increase in the extent of disturbance, and otherwise 
the solutions are almost identical. 

In the present section, the thin-layer solutions are 
compared with the solutions of the full Navier-Stokes equa- 
tions. In the first instance, at Re = 170, the difference bet- 
ween the thin-layer and full Navier- Stokes solutions is so 
small that it cannot be shown graphically. The profiles corr- 
esponding to this Reynolds number are already drawn in 
Figure 3.1* 

In Figures 5.15 and 5.16, the solutions of the thin- 
layer and full Navier-Stokes equations are compared at Reynol- 
ds numbers Re - 50 and Re = 10 respectively. As in the case 
of a cold wall, we observe that the thin-layer assumption 
gives thicker shock layers. The difference between the thin- 
layer and full Navier-Stokes equations increases predominantly 
as the Reynolds number is decreased. This is more, clearly 
brought out in Figures 5.17 and 5.18. In Figure 5.17, we draw 
the total extent of disturbance for the thin-layer and full 
Navier-Stokes equations, both for the slip case. Below Re -30, 
we observe that the two curves deviate enormously. At Re = 30, 
the thin-layer equations predict about a 10 percent thicker 
layer compared to the full Navier-Stokes equations. At Re = 20 
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thin is about 25 percent. At Re = 10, the thickness predicted 
by the thin-layer equations is almost twice the thickness 
predicted by the full Navier-Stokes equations. 

In Figure 5.18, we show the variation of the wall 
pressure with Reynolds number. Curves are drawn for the full 
Navier-S tokos equations and thin-layer equations, with slip 
as wo II as no-slip boundary conditions. The behaviour of the 
wall pressure for those cases is similar to that on a cold 
body, as depicted in Figure 4.21. The pressure prediction of 
the th.in-lf.iycr equations is not even qualitatively correct 
in the low Reynolds number regime. They predict a continuous 
drop of pressure* and do not show the rise at low Reynolds 
numbers. The full Navi or- St okes equations predict a steep 
rise at low Reynolds numbers, as observed in experiments 
(Potter and Bailey, 1963). The two curves for the full Navier- 
Stokes equations can be appended to the Section 5.3. Slip 
decreases the wall pressure. The effect of slip In the present 
insulated-wall case is more than that for a cold-wall case 
(cf 4 Figure 4.21 ). 

In Figure 5.19, we present the variation of slip 
velocity and the adiabatic wall temperature, with Reynolds 
number, for the thin-layer and full Navier-Stokes equations* 

In the thin-layer case, the wall temperature is always the 
same and is equal to the stagnation temperature of the free 
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stream. Thus the solutions of the thin-layer equations do 
not show the rise in the adiabatic wall temperatures at low 
Reynolds numbers as observed in experiments (for example, 

Drake and Backer, 1952). The full Navier-Stokes equations 
predict this rise. Regarding the slip velocity, Figure 5.19 
shows that the thin-layer equations over-est ima te it in the 
low- Roynolds-numb er regime. However, for Reynolds numbers 
greater than 30, the sslip velocity predicted by the thin-layer 
equations is slightly less than that predicted by the full 
Navior-S fcokee equations. This difference tends to zero as 
the Reynolds number increases. 

The above discussions lead to the following concl- 
usions regarding the validity of the thin-layer assumption 
with respect to the present adiabatic-wall solutions: 

(1) At low Reynolds numbers the thin-layer assumption 
gives solutions which deviate very much from the solutions 
obtained from the* full Navier-Stokes equations. 

(2) At low Reynolds numbers the thin-layer equations 
do not predict oven qualitatively the behaviour of certain 
quantities like the impact pressure and the adiabatic wall 

temperature. 

(3) At low Reynolds numbers the thin-layer equations 
over-predict the slip velocity and the extent of disturbance. 
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(4 ) At R ( -' — 50, the thin- lay er solutions show a maxi- 
mum difference of about 5 percent compared to the exact 
solutions. This difference increases fastly as the Reynolds 
number decreases, 

9 . 5 C omparison With Second-Order Boundary- layer Theory 

In this section we compare our results with the 
second -order boundary-layer theory for the stagnation point 
flow, Thu present results of the second-order boundary-layer 
theory are calculated from the solutions of Lenard (1962), 
who has integrated and tabulated the first- and the second- 
order boundary -layer equations, for various wall temperatures, 
for the stagnation-point flows around two-dimensional and 
nxisymuetrie bodies * 

In these comparisons, all our results presented 
are obtained from the solutions of the full Navier-Stokes 
equations lor an insulated wall. Further, we taJcey = 5/3, 

=-• 10 and Pr = ,75. A square-root viscosity-temperature 
law io used. The first-and second-order boundary-layer results 
ere calculated for corresponding parameters. 

In the first place, it may be noted that, according 
to the theory of lenard, the adiabatic wall condition 

(^I-) = 0 is satisfied, when and only when the wall temperature 

dr wall 
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is equal to the free stream stagnation temperature* This is 
also the case with the two -thin-layer theory of Cheng, as 
seen in Chapter 4 (cf. Figire 4.32). Thus both the second- 
order boundary-layer and the two-thin-layer theories fail 
to predict the experimentally observed rise of the adiabatic 
recovery temperature above the free stream stagnation tempe- 
rature, at low Reynolds numbers. However, as shown in 
Figure 5.7, the present full equations show this rise at 
low Reynolds numbers. 

In Figure 5.20* we present the skin- friction coeffl 


ciont comparisons. The results are given in terms of 


T Be oo 

t fm 



versus Re^,. The curves drawn are for the full Navier-Stokes 
equations with and without the slip condition, the first- 
order boundary- lay or theory, and the second-order boundary- 
layer theory, with all the second-order corrections included, 
as well as neglecting only the slip correction. Also the 
frce-molcculc limit for unit accommodation coefficients and 
infinite Mach number is shown in the figure. 

Fox* the present case, the second— order boundary— 


ley or theory gives the 
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T Rg oo 

T rcf 


1 ,3119 + 1 

ns* 

1st order boun- 
dary. layer 
theory 


curva- displa- slip vorticity 
ture cement 

- 1.603 + 2.07 - 2.927 + 1 ,672 j 

corrections due to 2nd-order boun- 
dary-layer theory. 


T ref iQ a reforGncG skin friction, which for the present 

case is related to the free-molecules skin friction 

2 

X FM ~ p °o V « sinO c °o9 'by 


T rof _ 3.762 



Thus wo get' 


t Re, 


other 

slip corrections 


oo 


4.937 _ 11.01 + 8.048 


M 


fRe 0 


Re_ 


Tie 


oo 


In the figure we have drawn two curves corresponding 
to the second-order boundary-layer theory; in one of them 
all the second-order corrections are included, in the other 
all the second-order corrections except the slip correction 

are included* 

Wo observe that, compared to the first-order boundary- 
layer theory, the second-order theory predicts a decrease 
of skin friction, when all the second-order corrections are 
included. But if the slip correction alone is neglected, then 
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the corresponding G^pve ri§§s above the fi^gt '-order boundary- 
layer curVe, fhe fl|ll Nqyier-Stokes equations algo show the 
saiiie trend regarding the effect of slip* When the slip is 
neglected the clirve is above the houndary-layer curve, and 
when slip conditions are Used the curve is below the bo unda ry— 
layer curve. But the magnitude of the slip effects differ 
considerably for the full Navier-Stokes equations and the 
second-order boundary-layer theory, particularly in low 
Reynolds numbers* Pot Re^ < 200:, the slip correction to skin 
friction provided by the full Navi er- Stokes equations is at- 
loast 60 percent larger than the slip correction provided by 
the second-order boundary-layer theory* However, at very high 
Reynolds numbers all the curves tend to values predicted by 
the first-order boundary-layer theory* 

In Figure we compare the slip velocity predic- 

ted by the full Navief^tokes equations with the secOnd- 
ordcr ho and ary-layer theory'* -For the present case the second- 
order boundary-layer theory gives the following expression 

■for the slip Velocity'. 



We observe from the figure that the agreement between the 
two curves is good at high Reynolds humbers^ Ehe slip Velocity 
from the second-order boundary-layer theory is always greater 
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than that predicted by the full Havier-Stokes equations. 

■The difiercnce between the "two curves gradually increases 
a.s the Reynolds number decreases, and below Re = 100 the 

.OO 

second— order boundary— layer curve shoots up very rapidly* 

In Figure 5*22, the tangential velocity profiles 
a_ k ^ G oo = ^00, 100 and 25 are drawn for our results and the 
second-order boundary-layer theory. It may be noted that 
accord ing to the second-order theory, the slip velocity on 
the body (shown in Figure 5.21') is contributed only by the 
slip correction. But in drawing the detailed profiles all the 
other corrections (curvature, displacement and vorticity) 
have their contribution. What we have drawn in the figure 
arc the boundary- layer profiles (Including the second-order 
corrections). This does not constitute the complete tangen- 
tial-velocity profile* The complete tangential -velocity pro- 
file, according to the second-order boundary-layer theory, 
should be calculated by correcting the inviscid flow due to 
the displacement thickness of the boundary layer, and combi- 
ning this Inviscid flow and the boundary layer profiles we 
have drawn in the figure, by some procedure of a composite 
expansion (Van Dyke, 1964a). However, very near the body sur- 
face the contribution of the inviscid flow is negligible and 
only the boundary layer solution dominates. 

For the "present values of the parameters, the boun— 
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a c o o I’d ing t o the Table 5,1, 

Table 5 . 1. Calculation of the Second-Order Boundary- 
layer Tangent ial-Velocity Profiles in 
Physical Coordinates. 


r 

first-order boundary layer + 

total second-order 




corrections to be 




applied. 

i 

0 

+ 

3.166 





.•JlliL. 

.4385 


2.591 

VKq*. 




1 .280 

. 6847 

+ 

1.804 





1.932 

.7881 

+ 

1.620 




T Re oo 


.8194 


.960 

Y»‘w 



f Re oo 


From Figure 5.22, 

we observe that at Re TO = 900, 

the prof, 

lies given by our 

integrations and the second-order 


boundary-layer theory are almost coincident upto some distance 
from the body. This distance constitutes about 15 percent 
of the distance from the body to the shock wave. But, at 
Re^ = 100, the 'two corresponding profiles slightly differ 
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from each other even on the body surface. And, at Re^ = 25, 
the difference increases sharply and the two profiles do 
not boar any relationship with each other. This shows the 
obvious failure of the second-order boundary-layer theory 
to predict the flow field at low Reynolds numbers. 

It is interesting to mention an anomalous behaviour 
of the HC'Oond-ordcr boundary-layer theory at low Reynolds 
numbers , when one attempts to obtain an uniformly valid solu- 
t.i on by cor hitting the second-order boundary-layer profile 
and the corresponding corrected inviscid profile. At Re^ = 50, 
it in found that the boundary layer (calculated upto the 
second order) happens to be so thick that it covers even the 
shock-wave position (duly corrected to account for the boun- 
dary layer thickness)* In such a situation it is impossible 
to obtain a meaningful composite solution. Similar anomaly 
is observed by Kao (1964b) with respect to his third-order 
boundary- lay or calculations. At Re^-60, he finds that the 
shock-wave thickness predicted by the third-order theory is 
so large that even a portion of the body is covered. These 
factors clearly show the failure of the higher-order boundary- 
layer theory at low Reynolds numbers. 

In brief, we can say that the second-order boundary- 
layer theory succeeds in predicting the qualitative behaviour 
of the overall characteristics of the flow such as skin 
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friction and slip velocity upto Re oo =30, even though the 
quantitative differences are appreciable. As far as the 
detailed structure ol the flow is concerned, predictions on 
the basis of the second-order boundary-layer theory are 
extremely poor in comparison to our results. The former even 
predicts physically implausible results below Re TO = 100. 

5 * 6 Comparison with Experiments 

In thin section, we present comparisons of our 
results with, the adiabatic -wall density profiles measured by 
Runsul (1968) and the impact pressure measurements of Sher- 
man (1993) and Potter and Bailey (1963). The density-profile 
comparisons for the cold-wall case have already been discu- 
ss,. .1 in Chr.pt or 4* In case of the present comparisons, Russel 
measured the density profiles with the body temperature equal 
to the f roe-stream stagnation temperature. As in the cold- 
wall cases, the f roe-stream stagnation temperature is kept 
equal to the room temperature (300°K). 
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Table 5.2: Test -Conditions for the Adiabatic -Wall 

Experimental Density Profiles With Which 
Comparisons are Made 


Ee m 

M 

oo 

Re 

s 

Re 

(Calculated from Re^ using linear 
viscosity-temperature law) 




Argon 

100 

3.83 

19.8 

16.97 

200 

3.83 

39.5 

33.94 

500 

3.8 

91.5 

85.97 

000 

3.8 

183.0 

171.94 




Nitrogen 

30 

4.2 

7.7 

6.62 

50 

4.38 

12.2 

10.34 

100 

4.38 

24.4 

20.67 

200 

4.19 

52.8 

44.34 

500 

4.11 

154.0 

144.2 


In Table 5.2, we give the conditions at which compa- 
risons are made. Both argon and nitrogen gas streams are 
used. The first three columns are taken from Russel's paper. 


205 


Here Re 


R=o Yoo 


s 


, where y corresponds to the viscosity 


calculated behind a normal shock wave* We have calculated 
Re from Re^ given by Russel by using a linear viscosity- 
temperature variation, and the corresponding values are tab- 
ulated in the fourth column. Some computations are also made 
using Sutherland viscosity law and it is checked that numeri- 
cal solutions in the two cases do not differ in any appreciable 


way. 


The comparisons are shown through Figures 5*23 to 5.31. 
The first four figures are for the argon profiles and the 
remaining for the nitrogen profiles. Numerical solutions, in 
all the cases, are obtained from the full equations with slip 
boundary conditions. The estimated accuracy for the experi- 
mental profiles is indicated in the figures. The profiles 
range from high to very low Reynolds number cases. The total 
extent of density disturbance ranges from about 25 percent 
to 65 percent of the body radius. The numerical solutions 
show a slightly thinner layer of disturbance than the experi- 
ments. From the figures, it is evident that the agreement 
is reasonably good. The general behaviour of the profiles 
in high as well as low Reynolds number cases is correctly 
predicted. 

Comparisons of our results with the impact-pressure 
measurements Sherman (1953) and Potter and Bailey (1963) is 
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presented in Figure 5.32, The present results are obtained 
integrating the full Navier-Stokes equatio ns with slip boun- 
dary conditions, andy = 1.4, =4,3 and 3^ = T^. The exp- 
erimental data of Sherman are for 1,8 < < 3,4 and of 

Potter and Bailey for 3.8 < < 5.8, with hemispherical 

probes at T w = T 0oo in a nitrogen stream. The impact pressure 
p^ is nondimensionalised with the inviscid impact pressure 
p^ which is related to p^, by 


p» 

o 

Pen 


y +1 


M 


y 

)Y r T 


2y 


Y+1 


M 


1/y —1 

jrnl) 


y+1 


The density ratio - 

, P s 1/2 

parameter ile s 


appearing in the correlation 


is calculated from the usual Rankin e- 


Hugoniot relation 


p p? _ y ~i 2 

p s Y+1 + (y +1 ) 55" 

The figure shows a good agreement of our results with the 
experimental data. As the Reynolds number is decreased the 
pressure first decreases from the inviscid limit and later 
steeply rises. The extensive measurements of Potter and 
Bailey (1963) for different flow conditions and probe 
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shapes show that this is the general behaviour of the impact 
pressures at lev/ Reynolds numbers. Our present results for 
various other body-temperatures and flow conditions conform 
to this behaviour (see Figures 4.21, 5.8 and 5.18). 

In Figure 5.33, the recovery factory, 

(I C)00 -Ico) , is compared with the experimental data of Hick- 
man. and Giedt (1963). The present results are obtained by 
solving the full Havier-Stokes equations with slip boundary 
conditions, fort= 5/3 and 1^=10. The experimental data are 
obtained using thin-walled glass models in an air stream at 
KU = 2-6. An important feature of the low-Reyno Ids -number 
results in" an increase of the adiabatic-wall temperature 
over the freo-stream stagnation temperature. This rise is 
due to the fact that in this Reynolds number range the vis- 
cous layer surrounding the sphere is so thick that heat is 
generated by dissipation more rapidly than it can be conduc- 
ted away. In spite of the differences in flow conditions, 
the present results show a reasonably good agreement with 
the experimental data. 
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. CHAPTER 6 

con erosions 

6,1 Introduction 

In this chapter, we epitomize the results arrived 
at in the present investigation. Most of these conclusions 
have already been mentioned in Chapters 4 and 5, where we 
have presented in detail the results of our work and compared 
them with several approximate continuum theories and avail- 
able experimental , results. Here, we shall reiterate the sali- 
ent features with a view to understand how far the questions 
raised in Section 1,5 are answered. These questions are: 

(1) What is the limit of validity of the If avier- Stokes 
equations in the hypersonic, rarefied flow around a 
blunt body? 

(2) How far does the thin-layer assumption give reasonably 
good results? and 

(3) What is the effect of the slip boundary conditions 
on the flow field? 

The present investigation is motivated by the fact 
that the continuum theories so far developed for the rare- 
fied blunt-body problem assume a thin layer of disturbance 
ahead of the body, and the Havier-Stokes equations are very 
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much simplified using this assumption. Solutions obtained 
from these approximate theories fail to predict several 
experimentally observed rarefied-flow phenomena. Jhirther*. 
all the existing continuum theories which include the slip 
boundary conditions, a priori assume a small effect on the 
flow field, and seek corrections to the no-slip solutions. 
Moreover, these corrections are calculated within the thin- 
layer model. In the present work, we have investigated the 
effect of the slip velocity and temperature jump using them 
as principal boundary conditions in the leading approxi- 
mation, and numerically integrated the resulting Uavier- 
Stokes equations near the stagnation region without making 
any further assumption regarding the model of the flow. 

We have obtained extensive numerical solutions for various 
degrees of rarefaction, wall temperature and gas properties. 
The important results are briefly discussed in the following 
sections, 

6.2 Continuum limit to the Hypersonic, Barefied 
Blunt-Bodv Problem 

It is long since known that the stagnation-point 
heat-transfer rates predicted by a simple and elegant 
continuum model, as in Cheng's two-thin-layer theory, . . 

agrees very well with the experimental values even in very 
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raroxied flow conditions. However, relatively recent 
experimental investigations (for example, Alt o use and Bog- 
donoff , 1969) on the detailed flow field in the stagnation 
region of a blunt body show that Cheng’s theory does not 
predict the detailed structure correctly even at Reynolds 
numbers as high as Re =60, Another recent approach to the 

o 

blunt-body problem is the use of direct Monte Carlo simu- 
lation, Vogenitz and Takata (1970) have compared their Monte 
Carlo profiles with the results obtained from the thin-layer 
continuum theory of Levinsky and Xoshihara (1962) at Re^ 152 
and reported poor agreement. The present continuum solutions, 
using the full Navier-Stokes equations, show a fairly good 
agreement with the above two cases, as shown in Sections 
4.7 and 4.8. The full Navier-Stokes equations predict the 
following additional rarefied gas-flow features that mary 
of the thin layer theories fail to give even qualitatively: 

(1) A reasonably good agreement of our theore- 
tical predictions with the available experimental, results 
for. the low Reynolds number flow profiles (Figures 4 #41 to 
4.49 and 5.24 to 5.31). 

(2) Agreement of our stagnation-region wall- 
pressure results with the experimentally observed behaviour 
(Figures 1.4, 5.8 and 5.32). 
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(3) Agreement of our theoretical predictions with 
the experimentally observed behaviour of the adiabatic wall, 
temperature in the low Reynolds number regime (Figure 5.33). 

These factors show that the full Ravier-Stokes equations 

(that is, without using any thin- layer assumption) reproduce 

many important features of the hypersonic, rarefied flow. 

The present investigation indicates that the Ravier-Stokes 

equations give physically plausible results even^a wall— 

Knuds en-number (Kr^ = ~) approximately equal to 0.1, or the 

B 

Reynolds number, Re, approximately equal to 3.0. Below Re=3.0, 
the heat-transfer results cross the free-molecule value and 
do not show any tendency to approach the free-molecule value 
as Re is further decreased. 

6.3 Thin-Layer Assumption at low Reynolds Numbers 

In the present study, the solutions of the full 
Navier-Stokes equations are compared with essentially three 
types of theories based on the thin-layer assumption. The 
first one is the numerical-integration-through-shock method 
applied to the thin-layer equations derived in Chapter 2. 

In the other two, namely, the hi^aen-order boundary-layer 
theory and the two- thin-layer theory of Cheng, the flow 
field is divided into distinct zones and is treated piece- 
meal. The main features of these three approaches are already 



described in Chapter 1 , where several arguments are put 
forward to show their inadequacy to describe very low 
Reynolds number flows . Retailed comparisons of the full 
Kav i er-S,t ok es equations with these approximate theories are 
made in Chapters 4 and 5. These comparisons conclusively 
show the failure of these theories at low Reynolds numbers, 
and throw more light on the nature of their failure. 

Comparisons of our full Navier-Stokes solutions with 
various approximate theories lead to the following conclusions 

(1 ) The thin-layer solutions differ' very much from 
our solutions in the detailed flow-structure* 

(2) The thin-layer solutions do not predict the 
correct behaviour of the stagnation-point pressures and the 
adiabatic-wall temperatures at low Reynolds numbers. 

(3) They predict lower maximum-temperatures in the 
flow field and predict far thicker merged layers than the 
full Navier-Stokes equations. 

(4) The extent of disturbance predicted by the thin- 
layer equations shoots up rapidly as the Reynolds number is 
decreased. The thickness of the merged layer at Re - 20 is 

of the order of the body radius, while the corresponding 
thickness from the full Navier-Stokes equations is much 
lesser. This .surprising situation, that at low Reynolds 
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numbers the solutions of the thin-layer equations predict 
thicker layers than the solutions of the full equations, is 
due to the cumulative effect of neglecting several viscous 
and heat-conduction terms in the equations. The neglected 
terms are comparatively unimportant at high Reynolds numbers, 
where there are large gradients in the viscous layer and 
the shock-wave-like regions. At low Reynolds numbers the 
corresponding gradients are small, and the neglected terms, 
which now become important, imply an underestimation of the 
viscous and heat conduction effects and hence predict thicker 
layers of disturbance. On the basis of overall comparisons, 
it may be said that the solutions of the thin- layer equations 
are not valid for Re < 50, having a difference of more than 
5 */. from the exact solutions. 

(5) A similar type of blowing-up of the profiles is 
observed in the higher-order boundary-layer theory as well* 
Here, it is found that at Re -20 the boundary layer becomes 
so thick that it even covers the shock-wave position, or the 
shock wave becomes so thick than it emcompasses even a portion 
of the body. Meaningful composite solutions cannot be obtained 
in such s ituat ions. 

(6) Both the higher-order boundary-layer theory and 
the two-thin-layer theory of Cheng do not predict the rise 
in the adiabatic-wall temperature and the stagnation-point 
pressure at low Reynolds numbers. 
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(7) Theories based on Gheng ! s model underestimate 
the effect of slip velocity and temperature jump. These 
theories give thinner layers compared to our solutions. In 
particular, the shock-transition zone is not correctly des- 
cribed by Cheng f s theory, further Cheng’s analytical results 
impose certain restrictions on the nature of the flow pro- 
perties like the use of a linear viscosity-temperature law. 

finally, we wish to make a remark on the use of local- 
similarity (Section 2.3). It has long been thought that the 
success of the local-similarity concept, to describe the 
stagnation-region flow, depends essentially on the concept 
of a thin layer of disturbance (e.g. Cheng, 1966, Mikhailov 
et al . , 1971). ^he extensive calculations and comparisons 
made in the present study, using local similarity in the 
full Ravier-Stokc-s equations and in the reduced, thin- layer 
equations, show the superiority of the former in describing 
several important features of the merged flows. This Indicates 
that the success of local similarity concept does. not intri- 
nsically depend upon the thin-layer approximation. 

6.4 Effect of Slip Boundary Conditions on the flow field 

(1) It is found that the slip conditions have any 
appreciable effect on the flow field only for Reynolds 
numbers, Re, less than about 300. Above this Reynolds number, 
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the slip velocity and temperature jump are very small and 
tend to zero as the Reynolds number is increased* 

(2) It low Reynolds numbers slip conditions at the 
surface influence the entire flow field including the shock- 
wave-like region. 

( 3 ) The combined effect of slip velocity and tempe- 
rature jump depends upon the wall temperature. Slip velocity 
and temperature jump have opposite effects on the flow field, 
Ror example, the slip velocity has a tendency to lessen the 
total extent of disturbance in the gas, while the temperature 
jump increases it. If the body is very cold, the effect of 
the temperature jump dominates over that of the slip velocity. 
At high wall-temperatures of the order of the free-stream 
stagnation temperature, the slip velocity effect dominates. 

The combined effect of the slip conditions on the flow field 
is more in an adiabatic-wall case than in a cold-wall case. 

(4) In a cold-wall case, the slip boundary conditions 

are found to increase the - skin friction and heat— transfer 

rate to the body, compared to the no-slip solutions. On the 

otherhand, in an adiabatic -wall case the slip boundary 

conditions decrease the skin friction. 

! ' 

(5) In both adiabatic and cold-wall cases, it is 

found that the slip conditions decrease the pressure at the 
stagnation point. 
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(6) It is observed that the effect of slip velocity 
and temperature jump intrinsically depends on the ratio of 
specific heats, y . For a larger y, the slip effects are 
larger. 

So far, no blunt-body experimental results are avail- 
able regarding the effect of the slip conditions, or the 
actual magnitudes of slip velocity and temperature jump. 

Hence comparison of our results with experiments is not 
possible. 

6 * 5 Some Problems for Further Research 

The success of the full Navier-Stokes formulation, in 
describing the merged-layer flow in the stagnation region 
of a blunt body, encourages the following important exten- 
sions using the present model: 

(1) Study of the nonequilibrium effects of dissocia- 
tion and ionisation. 

(2) Study of the radiative-transfer effects in the 
hypersonic, merged-layer flow about axisymmetric blunt bodies. 

(3) Investigation of flows with large injection 
from the body. 

(4) Investigation of the effect of higher-order 
series approximation about the stagnation line* 
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and 

(5) The present formulation is limited only to the 
stagnation region of a "blunt "body. These solutions may 
provide initial values for integrating the Havier-Stokes 
equations downstream of the stagnation point. 



1PPMDIX 


Table A-1 : Some Surface Characteristics of the Solutions 
of the Tull Favier-Stokes Equations with Slip 
Boundary Conditions^ Y = 5/3, 1^=10, Pr=3/4 and 
T w =0.029 


Re 

Kn 

w 

u(l) 

5(0 

TRe oo 

(v ^ 

Vic dr ; w 

1. p 

^OO oo 

5(i) 

t M 

2.55 

0.0799 

0.283 

0.452 

13.12 

1.186 

1.104 

3.40 

0.0691 

0.277 

0.438 

16.20 

0.988 

1.058 

4.25 

0.0568 

0.260 

0.414 

18.62 

0.821 

1.025 

5.10 

0.0416 

0.234 

0.363 

20.50 

0.809 

0.996 

6.80 

0.0247 

0.184 

0.293 

24.80 

0.691 

0.965 

20.40 

0.0046 

0.086 

0.154 

42.95 

0.428 

0.879 

25.50 

0.0033 

0.070 

0.131 

46.08 

0.371 

0.872 

34.00 

0.0019 

0.051 

0.105 

50.47 

0.307 

0.875 

42.50 

0.0013 

0.041 

0.092 

54.72 

0.269 

0.881 

51.00 

0.0010 

0.034 

0.083 

58.81 

0.242 

0.886 

59.50 

0.0008 

0.030 

0.075 

62.04 

0.220 

0.891 

68.00 

0.0006 

0*027 

0.070 

65.43 

0.202 

0.893 

76.50 

0.0005 

0.025 

0.067 

67.17 

0.185 

0.896 

85.00 

0.0005 

0.022 

0.062 

69.49 

0.173 

0.899 

102.00 

0.0004 

0.018 

0.058 

72.22 

0.152 

0.903 

136.00 

0.0002 

0.014 

0.052 

78.53 

0.123 

0.908 

170.00 

0.0002 

0.012 

0.047 

85.52 

0.107 

0.912 

221.00 

0.001 

0.009 

0.044 

88.15 

0.085 

0.914 

272.00 

0.001 

0.007 

0.040 

92.37 

0.072 

0.915 

296.00 

0.001 

0.007 

0.040 

98.87 

0.067 

0.916 
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Table A- 2 ; Some Surface Characteristics of the Solutions of 
the Pull Navier-Stokes Equations with No -Slip 
Boundary Conditions; Y= 5/3, = 10, Pr = 3/4 

and T w = 0.029. 


Re 


tR 6co 

cEr'w 

P (1 ) 

t pm 

1 p 

2 OO CO 

1.70 

0.0076 

7.89 

1.027 

1.280 

2.55 

0.0060 

9.35 

0.805 

1.178 

3.40 

0. 0043 

11.19 

0.728 

1.130 

4.25 

0.0036 

12.35 

0.640 

1.077 

5.10 

0.0030 

14.67 

0.604 

1.052 

10.20 

0.0016 

20.34 

0.429 

0*962 

13.60 

0.0012 

25.10 

0.394 

0.920 

17.00 

0.0010 

27.22 

0.344 

0.903 

20.40 

0.0008 

28.80 

0.305 

0.895 

25.50 

0.0007 

32.05 

0.270 

0.891 

34.00 

0.0005 

35.20 

0.224 

0.882 

42.50 

0.0004 

38.35 

0.193 

0.880 

51 .00 

0.0003 

43 . 60 

0.184 

0.894 

68.00 

0.0002 

50.17 

0.159 

0.897 

85.00 

0.0002 

54.93 

0.139 

0.899 

102.00 

0.00 02 

58.20 

0.123 

0.903 

119.00 

0.0001 

62.70 

0.113 

0.904 

170.00 

0.0001 

73.24 

0.092 

0.908 

187.00 

0.0001 

75.17 

0.086 

0.912 

238.00 

0.0001 

80.50 

0.073 

0.913 

255.00 

0.0001 

86.57 

0.072 

0.915 
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Table A-3 : Some Surface Characteristics of the Solutions 
of the Pull Navier- Stokes Equations with 
Slip Boundary Conditions ; y = 1.33? = 10, 



Pr 

= 0.71 

and T w = 

0.1. 



Re 


u(l ) 

T(1 ) 

xRe °o (k 

t fm T 

2 

0 

. dT% 

• 3r 'w 

TV 

oo v oo 

P(1) 

3.58 

0.0363 

0.187 

0.342 

12.64 

1.029 

1.199 

7.17 

0.0189 

0.152 

0.297 

18.06 

0,734 

1.062 

11.95 

0.0097 

0.113 

0.246 

26.21 

0.631 

1.025 

16.70 

0.0069 

0.096 

0.229 

30.99 

0.545 

0.982 

19.12 

0.0058 

0.089 

0.220 

33.15 

0.518 

0.965 

23.90 

0.0041 

0.079 

0.208 

35.89 

0.458 

0.952 

28.68 

0.0030 

0.065 

0.194 

39.40 

0.422 

0.949 

LT\ 

00 

• 

in 

0.0023 

0.058 

0.182 

42.26 

0.369 

0.947 

71.70 

0.0008 

0.035 

0.151 

56.43 

0.255 

0.929 

91.60 

0.0006 

0.029 

0.143 

62.86 

0.216 

0.927 

143.40 

0.0003 

0.022 

0.133 

71.96 

0.167 

0.911 

191.20 

0.0002 

0.017 

0.125 

79.10 

0.140 

0.915 
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Table 1-4; Some Surface Characteristics of the Solutions 

of the Pull Navier-Stokes Equations with No-Slip 
Boundary Conditions ; y = 1.33, = 10, Pr = 0.71 
and T w = 0.1. 


Re 


TR ®oo 

& ® ) 

^ dr'w 

P(1) 

t M 

2 ^00 ^00 

2.39 

0.0130 

9.92 

1.195 

1.469 

4.78 

0.0092 

13.73 

0.864 

1.254 

7.17 

0.0058 

17.87 

0.746 

1.090 

9.56 

0.0046 

21.01 

0.663 

1.040 

14.34 

0.0032 

25.43 

0.547 

0.987 

35.85 

0.0014 

39.49 

0.355 

0.939 

71.70 

0.0007 

52.28. 

0.242 

0.912 

119.50 

0.0004 

64.23 

0.181 

0.901 

143.40 

0.0003 

69.59 

0.164 

0.925 

239.00 

0.0002 

83.30 

0.120 

0.985 


Table A-5 : Some Surface 
of the Thin-] 
Condit ions ; y 


Re 

u(l) 

hi) 

11.9 

0.145 

0.198 

13.6 

0.123 

' 0.177 

17.0 

0.101 

0.156 

25.5 

0.068 

0.120 

34.0 

0.051 

0.099 

42.5 

0.040 

0.084 

51.0 

0.034 

0.076 

68.0 

0.027 

0.068 

85.0 

0.022 

0.061 

110.0 

0.012 

0.046 

255.0 

0.008 

0.042 
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Characteristics of the Solutions 
'er Equations with Slip Boundary 
5/3, M =10, Pr=3/4 and T =0.029. 

oo ■ W 


TRe 

OO 

(k & 

Vlc dr'w 

13 ( 1 ^ 

x m 

1 P v 5 

2 oo oo 

y v U 

24.34 

0.334 

0.605 

26.59 

0.332 

0.651 

31.55 

0.329 

0.706 

40.17 

0.292 

0.763 

45.28 

0.255 

0.789 

47.77 

0.229 

0.825 

52.40 

0.204 

0.830 

61.74 

0.174 

0.864 

66.40 

0.158 

0.887 

84.87 

0.104 

0.897 

96.70 

0.079 

0.912 


256 


Table A-6 : Some Surface Characteristics of the Solutions 

of the Thin- layer Equations with No-Slip Boundary 
Conditions ; Y = 5/3, M to = 10, Pr = 3/4 and T w = 0.029 


Re 

xR - e 0O 

ft S.) 

^ dr'w 

P(1) 

(1) 

t em 

6.8 

10.44 

0.192 

0.316 

52.73 

8.5 

12.41 

0.227 

0.457 

76.30 

10,2 

15.51 

0.265 

0.716 

119.44 

13.6 

20.24 

0.277 

0.793 

132,28 

17.0 

22.85 

0.256 

0.795 

132.54 

25.5 

28.05 

0.216 

0.795 

132.60 

34.0 

34.11 

0.200 

0.819 

136.68 

68.0 

47.04 

0.141 

0.843 

140.70 

85.0 

54.16 

0.131 

0.860 

142.00 

170.0 

76.13 

0.115 

0.920 

153.33 



Table it-7 : Some Surface Characteristics of the Adiabatic- 
Wall Solutions of the Pull Navier-Stokes 
Equations with Slip Boundary Conditions ; Y = 5/3 » 
= 10 and Pr = 3/4. 


Re 


50 ) 

TRe co 

t em 

T 

CO 

P(1) 

2.55 

0.381 

0.558 

13.59 

59.95 

1.017 

3.40 

0.297 

0.533 

15.69 

54.14 

0.883 

4.25 

0.281 

0.489 

18.38 

49.42 

0.844 

5.10 

0.192 

0.475 

20.06 

48.13 

0.815 

8.50 

0.097 

0.390 

24.55 

42.75 

0.834 

10.20 

0.082 

0.354 

27.97 

41.15 

0.843 

11.90 

0.077 

0.328 

31.07 

39.94 

0.846 

13.60 

0.063 

0.295 

37.92 

38.44 

0.850 

17.00 

0.043 

0.279 

41.94 

37.89 

0.855 

34.00 

0.020 

0.199 

59.20 

36.10 

0.870 

o 

o 

• 

in 

0.013 

0.171 

73.37 

35.77 

0.875 

85.00 

0.008 

0.128 

102.00 

35.08 

0,881 

93.50 

0.007 

0.123 

105.00 

35.04 

0.885 

102.00 

0.006 

0.117 

110.80 

35.02 

0.885 

119.00 

0.006 

0.110 

115.60 

35.00 

0.885 

136.00 

0.005 

0.105 

122.80 

34.92 

0.885 

153.00 

0.004 

0.096 

136.00 

34.80 

0.885 

170.00 

0.004 

0.093 

141.00 

34.70 

0.885 

187.00 

0.004 

0.085 

152.40 

34.57 

0.885 

221 .00 

0.003 

0.080 

163.15 

34.30 

0.885 

225.00 

0.003 

0.077 

167.50 

33.90 

0.885 
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Table A-8 : Some Surface Characteristics of the Adiahatic- 

Wall Solutions of the Full Navier-Stokes Equations 
With No-Slip Boundary Conditions; y = 5/3, 1^=10 
and Pr = 3/4. 


Re 

Kriyy 

v oo 

w 

3T 

oo 

poo 

t pm 

3.40 

0.198 

28.20 

48.36 

1.185 

4.25 

0,157 

30.72 

46.08 

1.119 

5.10 

0.138 

32.10 

45.12 

1.030 

6.80 

0.103 

36.00 

41.97 

0.983 

8.50 

0.085 

39.60 

41.03 

0.951 

10.20 

0.071 

42.50 

40.69 

0.926 

11.90 

0.060 

45.90 

39.46 

0.921 

13. CO 

0.051 

48.39 

38.44 

0.912 

17.00 

0.040 

53.08 

38.03 

0.910 

34.00 

0.019 

72.72 

36.58 

0.907 

51.00 

0.011 

87.47 

35.24 

0.902 

85.00 

0.008 

110.96 

35.21 

0.898 

93.50 

0.007 

115.57 

35.16 

0.890 

1 02 . 00 

0.006 

119.70 

35.04 

0;885 

119.00 

0.006 

129.58 

34.95 

0*885 

136.00 

0.005 

136.80 

34.92 

0.885 ' 

153.00 

0.004 

145.80 

34.75 

0.880 

170,00 

0.004 

152.10 

34.74 

0.880 

187.00 

0.004 

159.60 

34.46 

0.880 

221.00 

0.003 

170.60 

34.16 

0.878 

255.00 

0.003 

182.70 

33.80 

0.875 
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Table A-9 : Variation of Some Plow Characteristics With 

Wall Temperature; Pull Nay ier- Stokes Equations 
With Slip Boundary Conditions; Re = 6.8, Y = 5/5, 
= 10 and Pr = 3/4. 



u(1) 

1(1) 

T 

x Maz 

o€>w 

1 p v 

2 “ oo 

CD 

0.1 

0.210 

0.374 

0.792 

0.674 

13.10 

0.2 

0.238 

0.475 

0.825 

0.609 

9.95 

0.3 

0.264 

0.580 

0.872 

0.525 

8.02 

0.4 

0.294 

0.663 

0.895 

0.498 

6.75 

0.5 

0.313 

0.744 

0.945 

0.425 

5.97 

0.6 

0.319 

0.815 

0.958 

0.379 

5.42 

0.8 

0.349 

0.945 

1.038 

0.292 

4.50 

1.0 

0.376 

1.073 

1.135 

0.165 

3.95 

1.2 

0.404 

1.208 

1.210 

0.008 

3.46 
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Table ij --10 j Variation of Some Flow Characteristics With. 

Mach Number ; Pull Navier-Stokes Equations With 
Slip Boundary Conditions; Re^ = 40 , y = 5/3, 




Pr = 3/4 and T w = 

0.029. 


' K* 

oo 

u(1 ) 

5(i) 

%ax 

^ dr^w Reoo// Pooy 2 

OO 


4 • 0 

0.151 

0.335 

0.884 

15.25 

0.472 

6.0 

0.158 

0.290 

0.830 

19.16 

0.560 

8.0 

0.173 

0.286 

0.791 

21.75 

0.625 

10.0 

0.184 

0.293 

0.764 

24.80 

0.691 

12.0 

0.208 

0.323 

0.754 

26,42 

0.750 

■ 1 4 ! 

14.0 

0.232 

0.355 

0.752 

28.27 

. ■ . ■ " 1 

0.806 

16.0 

0.260 

0.395 

0.750 

29.54 

0.862 

18.0 

0.272 

0.415 

0.747 

30.24 

0.920 

20.0 

0.274 

0.424 

0.744 

31.08 

0.978 

22.0 

0.276 

0.429 

0.741 

32.10 

1.006 
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